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ABSTRACT 


DESIGN OF MULTIVARIABLE FEEDBACK CONTROL SYSTEMS VIA SPECTRAL 
ASSIGNMENT USING REDUCED-ORDER MODELS AND REDUCED-ORDER OBSERVERS 

Preston Ivanhoe Carraway III 

Roland R. Mlelke, Principal Investigator and 
Leonard J. Tung, Co-Principal Investigator 

The feasibility of using reduced-order models and reduced-order 
observers with eigenvalue/eigenvector assignment procedures is Investi- 
gated. A review of spectral assignment synthesis procedures Is present- 
ed. Then, a reduced-order model which retains essential system charac- 
teristics is formulated. A constant state feedback matrix vrfilch assigns 
desired closed loop eigenvalues and approximates specified closed loop 
eigenvectors Is calculated for the reduced-order model. It Is shown 
that the eigenvalue and eigenvector assignments made In the reduced- 
order system are retained when the feedback matrix Is Implemented about 
the full order system. In addition, those modes and associated eigen- 
vectors which are not Included In the reduced-order model remain un- 
changed In the closed loop full-order system. The full state feedback 
design Is then Implemented by using a reduced-order observer. It Is 
shown that the eigenvalue and eigenvector assignments of the closed loop 
full -order system remain unchanged when a reduced-order observer Is 
used. The design procedure Is Illustrated by an actual design problem. 
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CHAPTER 1 


INTRODUCTION 

The use of reduced-order models [1] and reduced-order observers [2] 
in the design of feedback controllers has been studied by several re- 
searchers. In addition, the development of eigenvalue/eigenvector as- 
signment techniques has received much attention in recent years. In 
this work, a reduced-order model is used with eigenvalue/eigenvector 
assignment techniques to design a constant state feedback controller for 
the original full-order system. The eigenvalues and eigenvectors con- 
tained in the reduced-order model are reassigned in the full-order sys- 
tem while those eigenvalues and associated eigenvectors not included in 
the reduced-order model remain unchanged in the full-order system. The 
constant state feedback matrix is implemented using output feedback with 
i a reduced-order observer. It is shown that the eigenvalues and eigen- 

vectors of the closed loop full-order system remain unchanged when the 
: reduced-order observer is implemented. 

1 . ' ' 

1.1 Motivation 

• During the past fifteen years significant advances have been made 

toward developing viable synthesis techniques for multivariable feedback 

i control systems. Notable among these techniques is the eigenvalue/ 
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eigenvector assignment procedure. Early studies In this area focused on 
an algorithmic formulation of the spectral assignment by Srlnathkumar 
[3], while later studies Included a geometric formulation of the same 
problem by Moore [4], Klmura [5], and Davison and Wang [6]. Based on 
these theories, design procedures have been developed for approximating 
desired mode mixing [7], reducing elgensystem sensitivity to variations 
In plant paroneters [8], reducing the effects of actuator noise on sys- 
tem performance [9] and modifying the resultant feedback gain matrix to 
specified gain constraints [10]. Recently, these procedures have been 
Incorporated In a spectral assignment computer aided design package 
[11]. A deficiency In all work concerning elgenvalue/elgenvector as- 
signment procedures Is an absence of application of these techniques to 
real world design problems. A primary factor contributing to this prob- 
lem Is the lack of understanding of how to use reduced-order models and 
reduced-order observers with spectral assignment procedures. 

Models representing the behavior of physical systems often consist 
of a very large nunber of coupled, linear differential equations. Such 
models are difficult to use when designing control systems due to exces- 
sive requirements for computer time and memory, and to the numerical 
analysis problems Inherently present when dealing with large systems of 
equations. It is, therefore, desirable to develop a design procedure 
which utilizes reduced-order system models. Simplification of large 
order dynamic systems h.’s received the attention of many researchers in 
recent years. The major difficulty with this work is that only open- 
loop system behavior is approximated. Of concern when using reduced- 



order models with eigenvalue/eigenvector assignment procedures Is the 
fact that while the reduced-order model may approximate open-loop system 
behavior, the modeling error may be so great or of such a nature that 
actual closed-loop system performance Is not acceptable. Also of con- 
cern Is the closed-loop behavior of those modes of the original system 
which are not Included In the reduced-order model . 

Full state feedback Is Implemented by the use of a dynamic observer 
system when there are fewer outputs than states. Since some states are 
usually available for measurement at the output, a reduced-order observ- 
er Is desirable In order to minimize the complexity of the control sys- 
tem. Of concern Is the effect of a reduced-order observer on the system 
eigenvalues and eigenvectors. 

1.2 Overview 

In this section an overview of the thesis Is given. A background 
of spectral assignment theory 1$ discussed In Chapter 2. A subsequent 
design procedure Implemented by Marefat In a computer aided design pack- 
age Is presented next. This Information provides a necessary foundation 
to support the material In the remaining chapters. In Chapter 3, a new 
technique Is developed that uses a reduced-order model of a known larger 
system and spectral assignment procedures to reassign selected eigen- 
values of the system, this Is accomplished without affecting the eigen- 
values and eigenvectors not included in the reduced-order model. Sec- 
ondly, a technique is developed that uses Luenberger* s [2] reduced-order 
observer and spectral assignment procedures to implement a constant full 



state feedback design using dynamic output feedback. The assigned 
eigenvalues and eigenvectors of the original system are retained using 
this technique. In Qiapter 4, a design philosophy and then a corre- 
sponding design procedure are developed for the new synthesis techniques 
presented In Chapter 3. A software package Is developed to facilitate 
the design of dynamic output feedback control systems using this new 
philosophy and procedure. The package Is Included as a new mode to a 
spectral assignment computer aided design program developed by Marefat 
[15]. The use of spectral assignment with reduced-order models and 
reduced-order observers In an actual design problem Is demonstrated In 
Chapter 5. Results are compared to those obtained by an alternate de- 
sign procedure. A program listing and an example of a computer aided 
design session are Included as appendices. 
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CHAPTER 2 

SPECTRAL ASSIGNMENT PROCEDURE 

In this chapter a background of spectral assignment theory Is pre- 
sented to support the development In Chapter 3. Definitions of eigen- 
values and eigenvectors are given. Then the effect of eigenvalues and 
eigenvectors on the time response of a system Is presented. Lastly, a 
characterization of the freedom available In selecting eigenvectors for 
a given eigenvalue assignment using constant state feedback Is pre- 
sented. 

2.1 System Elgenstructure and Time Response 
The eigenvalues of an nth order real matrix A are the zeros of the 
polynomial det [xi-A]. The eigenvalues, x^eC, form a self-conjugate 
set. That Is, for each complex eigenvalue Xf there exists a complex 
conjugate eigenvalue x <+ ^ » x^*. For each eigenvalue x^, there Is a 
right eigenvector, v^cC n , that satisfies the equation, 

A, f ■ V (2.1) 

for 1 ■ 1,..., n. If the eigenvalues of A form a distinct set, then 
the associated eigenvectors are linearly Independent [11]. Equation 
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(2.1) Is written for all A^ and v^ as 

AV ■ VA, (2.2) 

where V ■ [Vx ... V n ] and a ■ dlag . (xi ,...,A n )> Since the columns of 
V are linearly Independent, V Is Invertible. Therefore, 

A - VAV’ 1 . (2.3) 

Similarly, for each eigenvalue A^ there Is a left eigenvector, 

UjeC 1 ) that satisfies the equation 

Uj T A ■ *jUj T (2.4) 

for 1 ■ 1,..., n. The left eigenvector equation Is written for all A ^ 
and U| as 

U T A - AU T (2.5) 

where U * [ui ... u n ]. For distinct eigenvalues, the left eigenvectors 
are also linearly Independent [12]. Hence, premultlplylng equation 
(2.5) by (U T ) yields 

A « (U T )‘ l AU T . (2.6) 

Substituting for A from equation (2.3) Into equation (2.6) yields 





•1 . T. "I T 

VAV » (U 1 ) AU 1 . 


(2.7) 


Premultiplying by U and postmultlplylng by V yields 


U T VA ■ AU T V. 


( 2 . 8 ) 


Since A is a diagonal matrix of distinct eigenvalues, equation (2.8) 
can only be satisfied If U T V Is a diagonal matrix. For convenience 
the eigenvectors are usually normalized so that U T V ■ I or 


.1 


(2.9) 


The effect of eigenvalues and eigenvectors on system time response 
Is now presented. Consider the linear time Invariant system In Figure 
2.1 represented by the system state equations 


and 


x ■ Ax + Bu 


y ■ Cx, 


( 2 . 10 ) 


( 2 . 11 ) 


where A, B, and C are the plant. Input, and output matrices respec- 
tively and x e R n , u c R m , and y c R p . The system time response 
Is determined by solving the differential equation (2.10). Let a change 
of coordinates be defined by 
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x ■ Vz. 


( 2 . 12 ) 


The transformed system Is 


z ■ v’ l AVz + V -I Bu 


(2.13) 


y ■ CVz. 


(2.14) 


Substituting from equation (2.3) Into equation (2.13) yields 


z ■ Az + V Bu. 


(2.15) 


The solution of equation (2.15) Is given by [13] 

z(t) * e At Zq + {V (t “ T) U T Bu(t) dt (2.16) 

where zo Is the Initial value of z(t) at t * 0. Substituting from 
equation (2.12) yields the time response 


x(t) ■ Ve At U T x 0 + V I* e 


e A ^ U^Bu(tJ dT. 


(2.17) 


The first term of equation (2.17) Is celled the zero iiinfut response and 


the second term Is called the zero state response. 
Expanding the zero Input response yields 
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x(t) - [vi v n ] 


’ x.t A 

e 1 ... 0 

• 


T' 

Ul 

• • 


• 

• • • 



• • 



x t 

0 ... . e " 


s T 

■ « 


m m 


X0 


X,t. T . A , V . . 

v n e 1 (ui x 0 ) + . . . + v ln e (u n ) x 0 


Xit X t 

v ,e , T . V, T. 

nl (uj x„) + . . . + v e (u n ) x 0 


(2.18) 


From equation (2.18) the 1th component of the statevector Is determined 
to be 


n X *t j 

x.(t) « z v.. e (u. x 0 ). 
1 1-1 1 


(2.19) 


Expanding equation (2.19) yields 


Xi(t) 


vu 

Xit . T . 
e 1 (Ul X 0 ) 

v ln 

• 

• 

m 

i. 

m m 

+ . . . + 

• 

• 

• 

V 

nn. 


fhe zero state response Is expanded next. 


x t T 

e (u n x 0 ) 


. ( 2 . 20 ) 


u(t) be a vector of unit step functions, u 0 , for computational ease. 



ORIGINAL rv.v=" fg 
OF POOR QUALITY 


This yields 


x(t) - V / e A(t ’ r) U T B u 0 dT 


(Ve At ) / e _AT dx (U T B u 0 ). 


Since A is a diagonal matrix, the integral term is written as 


/ e" AT d. 


t , 

/ e XlT dx 


0 . . . . 


-XT 

l e n dt 


- -1 (e" Xlt -l) 


. . 0 


1 (e" Xnt -l) 


Premultiplying equation (2.22) by the diagonal matrix e yields 


hit -At 


-i (l-e Xl t ) 


... 0 


0 . . 


1 (l-e Xnt ) 
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The vector K Is defined to be 

K » U T B Uq . (2.24) 





» 

► 


s 


Substituting equations (2.23) and (2.24) Into (2.21) gives 


x(t) »[v! . . . v n ] 


-I (l-e^*) 


0 . . . 


... 0 


1 X n t 

- (1-e " ) 


n -K. x.t 
z Vi (-1) (1-e 1 ) 
1»1 1 x i 


(2.25) 


where denotes the ith element of K. Expanding equation (2.25) 
yields 


\(t) 


V 11 

(dk) (l-e^) 

v ln 

• 

> (t » 


m m 

Xl + + 

• • • • ' 

• 

• 

• 

V 

L nnj 


-K X t 

(— ) (1-e n ) 
Xn 


(2.26) 


x.t 

The terms e 1 are called the modes of the system. Equations (2.20) 
and (2.26) show that the eigenvalues of the system determine the rates 
of decay of the modes while the eigenvectors determine the contribution 
of each mode to the various states. Thus, the time response of a system 



can be controlled by proper selection of system eigenvalues and eigen- 
vectors. 

2.2 Characterization of Freedom in Eigenvector Assignment 
Given the linear, time invariant controllable system with constant 
state feedback in Figure 2.2, the system state equations are written 

x a (A + BF) x + Bv (2.27) 

y * Cx. (2.28) 

Given that constant state feedback is used, Wonham [14] states that an 
mxn matrix F can be found to assign an arbitrary self conjugate set 
of eigenvalues if the system is controllable. Moore [4] characterizes 
the freedom available to assign eigenvectors for an arbitrary self-con- 
jugate set of eigenvalues. He gives necessary and sufficient conditions 
to find a unique real matrix F that satisfies the eigenvector equa- 
tion 


(A ♦ BF) v, = V, A, (2.2,) 

For i = 1, ...» n when B has full column rank. Associate with each 
eigenvalue x^ an nx(n+m) matrix where 

i 
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<1 


and a compatibly partitioned (n+m)xn matrix 



s 



(2.31) 


r . 

t 


R 

I 


I 


t 


whose columns constitute a basis for the null space of S . Then the 

A i 

necessary and sufficient conditions to find a unique real matrix F 
that satisfies equation (2.29) are: 

1) Vectors v. eC n are linearly independent, 

. * * 

2) v- = Vj whenever X. * x. , and 

3) v e span (N, ). 

i A i 

Thus it is possible to assign an arbitrary selfconjugate set of 

eigenvalues and a set of eigenvectors from within the span of N. . The 

A i 

null space N. is determined by the selection of an eigenvalue X.. 

A 1 1 

The subspace, 1^ identifies the freedom available to assign eigen- 
vector V.j. 

2.3 Eigenvector Assignment for Real Eigenvalues 
It is first assumed that X. eR so that v^ e R n for i = 1,... n. 
Equation (2.29) is rewritten as 


(X.I-A) v.-(BF) v. 


= 0 . 



i 


t 


rz 


(2.32) 



WPWMl 


y' / 



i. 
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Since Kj^ Is a basis for the null space of S x , then any vector Kj 

that postmultlpl les 1C gives a resulting vector that lies In the null 

i 

space of S. . Therefore, 

A 1 


C^I-A : B] 


Expanding equation (2.33) yields 


"a, 

«Al 


K i - 0. 


(2.33) 


[A^-A] N^K i + [B] M xi K i - 0. 


(2.34) 


Since Vj e span (N^), then determines where In the allowable 


subspace v^ exists. Hence, 


v i - % K r 


It follows from equations (2.32) and (2.35) that 


F. ( - - K,, 


By defining w^ as 


(2.35) 


(2.36) 


equation (2.36) Is rewritten in matrix form for all 1 as 


(2.37) 
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F[vj t •••• v^] ■ [w lf (2.38) 

or 

FV*W. (2.39) 

Since the eigenvectors are linearly Independent, then 

F « W V- 1 . (2.40) 

2.4 Eigenvector Assignment for Complex Eigenvalues 
It Is next assumed that A^eC for 1 ■ 1, 2 and A^eR for 
1 »3, ..., n. Then the first closed loop right eigenvector must satisfy 
the equation 

[A + BF] (v RE +jv I|y| ) * (v RE + jv IM ) (A re + JA im ) 

t 

where the subscript one Is suppressed for simplicity. Equating 
Imaginary parts yields 

[A ♦ BF] v R£ ■ v R£ X RE - v IK X, H (2.42) 


(2.41) 
real and 


and 



18 


[A + BF] v IM ■ v IM X RE + v RE X IM . (2.43) 


The two equations are written In matrix form as 


and 


[X re I-A ! X IM I i B] 


RE 

■ V IM 

' Fv rei 


(2.44) 


Ore 1 '" 


x im ! 


B) 


'IM 

hV RE 

-Fv 


IM 


0 . 


(2.45) 


For the case of complex eigenvalues, the nx(2n+m) matrix Is defined 
as 

S X * O re I-A:X im I:B] (2.46) 

and a compatibly partitioned (2n+m)xn matrix *K^ Is defined by 



(2.47) 


where the columns of constitute a basis for the null space of Sy 
Hence, 
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Cx rb i ‘ a!a im iJb] 


(2.48) 


From equations (2.44), (2.45) and (2.48) It Is apparent that the vectors 
In (2.44) and (2.45) are contained In the null space defined by K^. 
Therefore 


and 


•m 

V RE 



" V IH 

e SPAN 

p x 

; Fv re 


M x 

■ m 


(2.49) 


V IM 


"x 

V RE 

e SPAN 

P X 

~ Fv im 


M X 

p • 


» « 


(2.50). 


From equations (2.49) and (2.50) It Is apparent that the allowable sub- 
space for v R £ and Vj^ Is described by 



e SPAN 



(2.51) 


and 


I 
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RE 

'lM 


'e SPAN 


(2.52) 


Combining equations (2.51) and (2.52) yields 


'RE 

'iM 


e SPAN 


-P, 


r\ SPAN 


N, 


(2.53) 


The following characterization of the freedom available to assign 
complex eigenvectors is not developed here but Is proved by Marefat 
[15]. Matrixes a and P are defined by 


-P. 


(2.54) 


and 


(2.55) 


K and K fl are defined to be matrices whose col inns constitute bases 

a p 

for the null spaces of a and p. respectively. Matrix y Is defined 
by 


T * CK 0 ‘- V 


(2.56) 


rz. 
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and 1^ Is defined to be a matrix whose columns constitute a basis for 
the null space of y. .A basis for [SPAN(a)ASPAN(e)] Is [(SPAN(o)) 1 + 
( SPAN(e )) J '3‘ L where denotes set direct summation and "i M 
denotes orthogonal complementation. Also, a basis for y 1 Isa basis 
for y T . Hence 

e SPAN (K^). (2.57) 

A specific vector within the null space of y Is defined by postmultl- 
plying by a vector x^. Thus 




K Y *V 


(2.58) 


Using equations (2.51) and (2.52), x c ^ and x^ are defined by 



(2.59) 


(2.60) 


The. left Inverses of 
linearly Independent. 



Therefore 


exist since the columns of 


are 
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From equations (2.48) and (2.49) It Is apparent that the vectors [-Fv RE ] 
and C-Fvjm] He In the space defined by the columns of M^. Hence 

- Ftf RE ■ \ x cl <2 - 63) 

and 



M A x c2* 


(2.64) 



Since Aj and \ 2 e C and *1 e R for 1 ■ 3, ...» n, then 
★ 

*2 ■ because the eigenvalues form a self conjugate set. Further- 
more, the second condition of spectral assignment requires that v 2 * 

* 

Vj . Thus the specification of one complex eigenvalue and eigenvector 

contains all the essential Information of the complex conjugate pair. 

★ 

It Is also important to note that If v 2 » v 2 and the pair v 1# v 2 are 
linearly Independent, then v RE and v^ M are also linearly Inde- 
pendent. In order to calculate the feedback matrix F, the following. 




1V2 

'J 

A ' 


definitions are given: 



Wl - M^, 

(2.65) 


W 2 - -M x X c2 , • 

(2.66) 


vi -v RE , 
and 

(2.67) 


V * V IH* 

Recalling that for the case of real eigenvalues 

(2.68) 


», • -M^ x, 
and 

(2.69) 


v i"V 

(2.70) 


equation (2.38) Is rewritten so that 



F ^ v RE’ v IM’ v 3* •••» v n^ ’ ^“ M \ x cl*“ M X x c2*" N X 3 x 3’ ,,,, ‘ M X n x n^ 

(2.71) 


Substituting equations (2.65) through (2.70) Into (2.71) yields 


/ 

F[v lt .... v n ] « [wj, .... w n ] 

(2.72) 

» 

or 



FV - W. 

(2.73) 


a 


i 
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As In the case for real eigenvalues, 

F » WV- 1 . (2.74) 

This development Is easily extended to more than one pair of complex 
conjugate eigenvalues. 

2.5 Use of Eigenvector Freedom 

It Is shown In the previous three sections that eigenvectors for 
the selected eigenvalues must reside In an allowable subspace that Is 
determined by the plant matrix A, the input matrix B, and the se- 
lected eigenvalues a^. Normally the eigenvector assignment that is 
most desirable for a given set of eigenvalues Is not achievable because 
It does not lie within the allowable eigenvector space. In this case It 
Is desirable to select the allowable eigenvector that is closest to the 
desired eigenvector. This Is accomplished by projecting the desired 
vector Into the allowable space so that the error between the desired 
and the assigned vector Is minimized In a least squares sense as Illus- 
trated in Figure 2.3. 

The desired vector is projected onto the allowable space by the 
projection operator 

<2 - 75) 

for real eigenvalues and 
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p Ky * k y OSf\)’V 


(2.76) 


for complex eigenvalues £15]. Indicating the desired vector by the 
subscript "D" and the assigned vector by the subscript "A", the 
projection is accomplished by the equations 


and 



V RE 

It 

V RE 

V IM 

A 

V !M_ 


(2.77) 

(2.78) 
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for the real and complex assignments, respectively. 

2.6 Improvement of Initial Assignment by Gradient Search 
The selection of eigenvalues and eigenvectors for a system is nor- 
mally motivated by the desire to shape the time response as discussed in 
Section 2.1. However, once the desired time response is approximated, 
there are often other aspects of the assignment that are unacceptable. 

An example is an assignment which requires extremely high feedback gains 
which are expensive to implement and very sensitive to noise. Another 
example is extreme eigensystem sensitivity to small plant parameter 
variations or modeling errors. The freedom available to select the 
eigenvectors often provides a means to drastically improve these seconJ- 
ary design objectives while only slightly modifying the initial eigen- ' 
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vector assignment and thus the time response. This improvement is ac- 
complished by modifying the eigenvectors within an area local to the 
original assignment. The vectors are modified in such a manner as to 
reduce the undesirable aspect of the assignment most rapidly. 

A cost function J is defined so that a reduction in the value of 
J corresponds to reduction of the undesirable aspect of an initial 
eigenvector assignment. A gradient matrix is computed in terms of J 
to determine how the eigenvector assignment is most efficiently changed. 
Recalling that the eigenvectors are determined by the equation 



for the case of real eigenvalues and 


v 1re 

* K 

i 

v ilM 

Y 



(2.79) 


(2.80) 


for the case of complex eigenvalues, it is apparent that small varia- 
tions in X. will cause correspondingly small variations in the eigen- 
vector assignment. A matrix X is defined as 

X - [X lt .... X n ]. (2.81) 


Since this matrix designates which eigenvectors are assigned, it is 
called the designator matrix. A gradient matrix [GR] with elements 
[GR]^ is defined to be 


l: 
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CGR]^ 


30 



M ~-H 

3X ij 


(2.82) 


The designator matrix X is then varied according to the rule 

X ij (q+1) * X ij (q) “ d[GR] ij (2.83) 

where d denotes the step size of X.. during each iteration. The 
gradient search is continued until a satisfactory compromise between the 
reduction in the value of the cost function and the modification of the 
time response is achieved. 



CHAPTER 3 


SPECTRAL ASSIGNMENT USING REDUCED-ORDER 
MODELS AND REDUCED-ORDER OBSERVERS 


In this chapter, the use of reduced-order models and reduced-order 
observers in the design of feedback controllers is investigated. A 
reduced-order model of a known system is formulated. It is then used to 
design a constant full state feedback matrix for the original full-order 
system. It is shown that the eigenvalues and eigenvectors reassigned in 
the reduced-order model are reassigned in the full-order system while 
those not included in the reduced-order model remain unchanged. The 
constant state feedback matrix is then implemented by output feedback 
using a Luenberger [2] reduced-order observer. It is shown that the 
eigenvalue and eigenvector assignments in the full-order system remain 
unchanged when a reduced-order observer is used. 


3.1 Motivation for Using Reduced-Order Models 
and Reduced-Order Observers 

Models representing the behavior of physical systems often consist 
of a very large nunber of coupled linear differential equations. Such 
models are difficult to use when designing control systems due to exces- 
sive requirements for computer time and memory, and to the nunerical 
analysis problems inherently present when dealing with large systems of 
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equations. It is, therefore, desirable to develop a design procedure 
which utilizes reduced-order system models. 

The spectral assignment synthesis methods described In Chapter 2 
use full state feedback. However, full state feedback Is not feasible 
for most systems because there are often fewer outputs than system 
states. Full state feedback is Implemented by the use of a dynamic ob- 
server for these systems. The use of a full system observer is unneces- 
sary since some states usually are available for measurement and there- 
fore need not be estimated. A reduced-order observer Is therefore de- 
sirable In order to minimize the complexity of the control system. 

• This chapter develops a reduced-order model and reduced-order ob- 
server. Control system design for the full-order system Is accomplished 
uslnq the reduced-order observer. Reduced-order models and observers 
have been used for several years. However, It Is shown here that the 
eigenvalues and eigenvectors assigned using the reduced-order model are 
retained In the closed loop full-order system while the eigenvalues and 
eigenvectors not Included In the reduced-order model remain unchanged in 
the closed loop full-order system. 

3.2 Reduced-Order Model Formulation 
A reduced-order model that Is used in the design of a constant 
state feedback controller for the full-order system model Is formulated 
in this section. The reduced-order model contains the eigenvalues that 
are to be reassigned In the full-order model. Let the original system 
model be described by the state equations 




and 


x - Ax ♦ Bu (3.1) 

y - Cx (3.2) 


where x eR n , ueR m , and >cR p . It Is assuned that the eigenvalues 
of A are distinct and are denoted by (Aj, ...,x n ). The corresponding 
modal matrix for A Is denoted by V * [v lt ...» v n ] where v^ de- 
notes the eigenvector corresponding to A ^ . The system model Is trans- 
formed by defining a new state variable 

z • V -1 x. (3.3) 


Equation (3.1) Is transformed to give 


z ■ Az ♦ Bu (3.4) 

-l , -l 

where A * V AV * diag (x 1# . .., x n ), and B » V B. The system Is 

now partitioned to separate the eigenvalues to be reassigned In the 

reduced-order model from those that will remain unchanged In the full- 

order system model. Thus, 



(3.5) 



j 


a 

4 


U 


It n a |r 

where z^R and Z 2 CR . The eigenvalues Included In the reduced- 
order model must be contained In Aj and those which are not Included 
must be contained In A 2 . The reduced-order model Is thus described by 
the state equation 

z l “ Ai Zi + Bi U. ( 3<6 ) 

The reduced-order model Is then used In conjunction with the spectral 
assignment procedure to assign eigenvalues and partially assign eigen- 
vectors in the full-order system model. However, the relationship 
between the eigenvalue and eigenvector assignments In the reduced-order 
and full-order models must be Investigated first. 

3.3 Spectral Assignment Using Reduced-Order Models 
The reduced-order model Is used to design a constant state feedback 
matrix for the full-order system. The relationship between the eigen- 
values and eigenvectors of the closed-loop reduced-order model and the 
closed loop full-order system must be understood In order to accomplish 
this. The relationship between reduced-order and full-order system 
eigenvalues Is determined first. Let F denote a constant state feed- 
back matrix computed for the reduced-order model. The control law Is 
then written as 


u « F zi. 


(3.7) 


The reduced-order model closed loop equation Is therefore 


A £J 


/ 

/* 


tn 

i 

i 



2 i • (A1 ♦ §1 F) Zi. ' (3.8) 


F Is now liripl emented about the full-order system by assuming full state 
availability in the full order model and transforming F back to the 
original coordinate system. Equation (3.7) is rewritten as. 


u ■ [F:0] “ F'z. 

2 2 

Substituting for z from equation (3.3) yields 

-i - 

u * F'V x « Fx. 

Hence, the closed-loop full order system is written as 


(3.9) 


(3.10) 


x * Ax + B Tx 
■ [A B F]x. 


Oili) 


The eigenvalues of the full order system are the eigenvalues of 
[A + BF]. This matrix is rewritten as 

[A + BF] « [VAV _l + VBF'V -1 ] 

» V[A + BF* 3 V" 1 . (3.12) 


« BBI UM 
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Since [A+BF] and [A + BF'] are related by a similarity transform- 
ation, they have the same eigenvalues. Matrix [A + BF 1 ] Is expanded 
as 


[A ♦ BF' ] ■ F-^-l + [Ii&qI 

[o:a 2 J [izF-oJ 


lAutJlIlfll 

I A A # I 

L 02F :A2 J 


( 3 . 13 ) 


The eigenvalues of this matrix are obviously the eigenvalues of 
[Ax + Bi F] and [a 2 ]. Thus It Is possible to reassign the k eigen- 
values included In the reduced-order model without modifying the (n-k) 
original system eigenvalues which were not Included In the reduced-order 
model. 

The relationship between eigenvectors of the reduced-order and full 
order system models Is determined next. The eigenvector equation for 
[A +BF'] Is written as 


L -b 2 f : *,-a 2 J V 2 , 


« 0 ( 3 . 14 ) 


which yields the two equations 


[A 1 I-Ai -Bi?} v lt » 0 


and 


I -M3 + [x 1 i-a 2 ] v 2i - 0. 


( 3 . 15 ) 

( 3 . 16 ) 
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If A^ Is an eigenvalue of [A t + §! F], but not an eigenvalue of 
[a 2 ], then from equation (3.15) It follows that v^ Is an eigenvector 
of [Aj + BjF]. Since \ ^ is not an eigenvalue of [a 2 ], [ajI-a 2 3 1s 
nonsingular. Thus, from equation (3.16) 

v 2i » (A i I-A 2 ]" 1 B 2 Fv 11 . (3.17) 

Therefore, v^ Is written as 


K 


- 




[a^i-a 2 ] 




(3.18) 


Equation (3.18) illustrates that the first k elements of v^ can be 
assigned using the reduced-order model while the remaining (k-p) ele- 
ments are linear combinations of v^ . 

On the other hand, if A^ Is an eigenvalue of A 2 and not an 
eigenvalue of [A x + BjF], the matrix [a^I-Ai-I^F] Is nonsingular. 
Therefore, equation (3.15) Is statlsfled only If 


Vjj ■ 0. .(3.19) 

From equation (3.16) It follows that v 2 ^ Is an eigenvector of A 2 for 
eigenvalue A^. Since A 2 Is a diagonal matrix, v 2 ^ Is written as 

$21 " £0* •••* 0* •••» (3.20) 
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where k^ Isa nonzero constant. 


Therefore v^ Is written as 



(3.21) 


Let Vj be the eigenvector of [A+ BF] corresponding to eigenvalue 
A j . The eigenvectors 7j and v^ are related by the transformation 



V v 


r 


(3.22) 


Expanding equation (3.22) for eigenvalues of [Ai + BjF] yields 


7 11 

9 

V u : v i 2 


V 21 


V 21 : V 22 

■ ■ 



V 11 + V 12 - 
V 21 V 22 CV’V’Vj 11 ' 


(3.23) 


Often when the reduced-order model contains only the dominant modes 
of the system, Vn * Vn + Vn [a^I-A 2 ] B 2 F. Then the top k compo- 
nents of the first k eigenvectors are assigned by choosing v^ as 


5 n - v n V 


(3.24) 
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When tire above approximation does not apply, an Initial assignment of 
0^ Is made using equation (3.24), F Is calculated, and the error 
between the desired top k components of the first k eigenvectors and 
the actual assignment Is calculated. A gradient search procedure Is 
then used to reduce this error. 

For Aj which are eigenvalues of A 2 , equation (3.22) expands to 

7 i * V5 1 ■ [&] ? 21 • (3 ’ 25) 

Substituting for from equation (3.20) yields 




(3.26) 


Therefore, the last (n-k) eigenvectors In the closed loop full order 
system are the original open-loop eigenvectors. Thus eigenvectors cor 
responding to the eigenvalues retained In the full order system model 
remain the same In the final closed loop system. 


3.4 Reduced-Order Observer Formulation 
In order to Implement the feedback matrix F calculated In equa- 
tion (3.10), full state availability is required. However, If the 
number of outputs p Is less than the number of states n, then (n-p) 
states must be estimated. This section parallels Luenberger's [2] de- 
velopment of a reduced-order observer system to allow the Implement a- 
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tlon of a full state feedback matrix In such a system. The following 

, section shows that not only are the eigenvalues of the observer system 
# 

retained In the closed loop system, but that the original eigenvector 
assignment is not affected by use of the observer. 

The open loop full-order system model described by equations (3.1) 
and (3.2) Is 

x ■ Ax + Bu (3.26) 

and 

y-Cx. (3.27) 

If the feedback matrix F is Implemented about the system model, then 

u-Fx. (3.28) 

It Is assumed without any loss In generality that the first p columns 
of C are linearly independent. A transformation matrix M Is defined 
to be 

■f u 

i.L 

where I n _ p Is the (n-p) order Identity matrix, 
is defined by 


(3.29) 

A new state variable 
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t 

z ■ Mx. 

(3.30) 

Equations (3.26), (3.27), and (3.28) are transformed 

and written as 

• 

M iy 

z * A z + Bu, 

(3.31) 

y-*?. 

(3.32) 

and 


«»p; 

(3.33) 

where 


A » MAM -1 , 

(3.34) 

B - MB, 

(3.35) 

l ■ CM' 1 • [I p : Oj, 

(3.36) 


and 
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~ - _i 

F » FM . 

Equations (3.31) and (3.32) are expanded as 


(3.37) 


mt 

• 

*1 


^11 

a 12 


Zl 


Bi 


3 





+ 


A 



1 

CM 

CM 

><C 


.*2. 


.V. 


(3.38) 


y * [Ci 




(3.39) 


Therefore the p components of the output vector are the first p 
states of the transformed system denoted by Zi . The reduced-order 
observer must then estimate the remaining (n-p) states denoted by Z 2 . 
Equation (3.38) is expanded to be 


and 


zi 


Z2 


An 

Zl 

CM 

+ 

Z2 

+ Bi u 

(3.40) 

A 21 

Zl 

+ A 22 

Z2 

+ B 2 u. 

(3.41) 


Since \ Is available as the output vector y, it can be different- 
iated to generate Zj . Hence, equation (3.40) is solved for A 12 z 2 
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which Is used as an input to a reduced-order observer to approximate 
z 2 . The proposed observer is shown in Figure 3.1. It is desired that 

e - z 2 (3.42) 

in order to implement the feedback matrix F. The observer state equa- 
tion is written as 

e = E6 + Lz t - LA n Zj + [n-LBju. (3.43) 

The need to differentiate zi is avoided by redrawing the observer 
as shown in Figure 3.2. If w is defined by 

w ■ 0-LZi , (3.44) 

then 

w 3 Lz - 6 3 0. (3.45) 

Using equation (3.42) to substitute for 9 in equation (3.45) yields 

• • 

w + L Z; - z 2 = 0. (3.46) 


Calculating each term of the above equation results in the three equa- 
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original; f*£ ?.v2 , 
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tlons, 

M ■ E [w ♦ Lzi 3 ♦ [fl-LBi]u - LAnZi, (3.47) 

III * LAj x + I-A 12 22 + LBi u, (3.48) 

# 

and 

-2 ■ -A 21 Zi - A 2 2 Z2 ' B 2 U. (3.49) 

Substituting equations (3.47), (3.48), and (3.49) Into (3.46) gives 

Ew ♦ [El-A 2l ] Zi ♦. [LAi 2 -A 22 ] z 2 ♦ [Q-B 2 ] u * 0. (3.50) 

Using equations (3.42) and (3.44) to substitute for w yields 

[A 21 ] Zj + [E-A 22 + LA 12 ] z 2 + [n~I 2 ] u * 0. (3.51) 

Since the Input and state vectors are not generally zero, the multiply- 
ing matrices must all be equal to zero for the equation to be true. 
Solving for the last two terms gives 

E ■ A 22 - LA 12 (3.52) 
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* 

and . 

ft«B 2 . (3.53) 

« 

Matrix Is generally nonzero also. This indicates that the pro- 
posed reduced-order observer is not adequate. Thus the observer Is 
modified by adding ’k zl Z x to w* and grouping terms as shown in Figure 
3.3. To remove the simmer located after the integrator, it is noted 
that 


EL z, » [(A 22 - LA 12 )L] Z} . (3.54) 

Using equation (3.54), the reduced-order observer is drawn as Figure 
3.4. The eigenvalues and eigenvectors of E ■ [A 22 - LA l2 ] are deter- 
mined by proper selection of L since the eigenvalues of 
^A 22 - LA 12 ] are also the eigenvalues of [A 22 *- A J2 T L T ]. Chapter 2 
describes a procedure for selecting a proper iJ to achieve a desired 
eigenvalue and eigenvector assignment. A guideline for selecting 
reduced-order observer eigenvalue locations is discussed in Chapter 4. 

The reduced-order observer is now used to implement the feedback 
matrix F. Expanding the control law given by equation (3.33) gives 

u * [Pj: ? 2 ] 


r~ ‘ 
Zl 


F 1 2 1 ♦ F 2 2 2* 


(3.55) 
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Using equations (3.42) and (3.44), a substitution Is made for z 2 re- 
sulting In 


u ■ CFi ♦ F 2 L] Hi * [F 2 ]w. (3.56) 

The control law In equation (3.56) Is Implemented In Figure 3.5. 

3.5 Effect of Reduced-Order Observer on Full System Elgenstructure 
The effect of a reduced-order observer on the eigenvalue and eigen- 
vector assignment In the closed loop full -order system Is developed In 
this section, luenberger [2] has proven that the eigenvalues of the 
original system assignment and the observer assignment remain unchanged 
in the closed loop full-order system. It is shown here that the eigen- 
vector assignment also remains unchanged. The following definitions for 
matrices 6, R, and T are given to reduce the algebraic complexity of 
this development. Let 


C(A 2 i - LAj i ) + (A 22 - LAi 2 )L], 

(3.57) 

R - [Fi + F 2 L] 

(3.58) 


and 


T - C-L: I]. 


(3.59) 
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The eigenvalues of the system are determined first. Using equa- 
tions (3.57), (3.58), and (3.59), the total system state equation Is 
written as 


A ♦ BR& BF 2 

TBRC + *GC E + TBF 2 


(3.60) 


Equation (3.60) Is now transformed to an upper triangular form so that 
the system eigenvalues are apparent. The transformation matrix P Is 
defined by 


[.: :]• 


and a new state variable Is defined to be 


(3.61) 


v ■ w- Tz 


(3.62) 


so that 




■ *1 " 

2 

s 

1 0 

• 


.-T ij 

M 

L v J 



BRC BF 2 

TM + 615 E + tBP s 


;] [;] • 


(3.63) 


Equation (3.63) Is simplified to give 


i 


0 


/ ' 
/• 
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/ / 


A + B (RC + F T) BF 2 
-TA + (£ ♦ ET E 
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(3.64) 


The submatrix [-TA + <& + ET] Is equal to zero. This Is shown by eval- 
uating each term within the matrix. The Individual terms are given by 


-TA ■ [LA,, - A 2 , • LA 12 - A 22 ], 


(3.65) 


GC - G[I:0] » [A 2 , - LA,, + (A 22 - LA 12 ) L : 0], (3.66) 


and 


Hence, 


ET . C-(A 22 - LA 12 )L : -LA, 2 + A 22 ]. 


(3.67) 


-TA ♦ {£ + ET « 0. 


(3.68) 


The expression (RC + F 2 T) Is equivalent to F. This Is shown by ex- 
panding (RC + F 2 T) as 

RC + F 2 T » [F, + F 2 L] [I : 0] + F 2 [-L : I] 


CFj : F 2 ] - F. 


(3.69) 




/ 


Therefore the substitution of equations (3.68) and (3.69) into (3.63J 
gives 


A + BF BF, Z 


(3.70) 


Thus the eigenvalues of the system are those assigned to £a + BF] and 

[E], In other words, the use of the reduced-order observer has no 
effect on the original system eigenvalue assignment. 

The eigenvectors of the system are determined next. Transforming 
the state equation using equation (3.30) results in 


* // 


m o huw m 


0 I 


8 ? Bl 


][:) 


M 0 |x 


01 I v 


(3.71) 


Let the closed loop system matrix in equation (3.70) be denoted by A r 
An eigenvector of corresponding to X Is denoted by Vp Eigen- 
vector v T is compatibly partitioned so that 


" n 







(3.72) 


The eigenvector equation for Ay is 


A I • Aj - BF 2 I fv 


*1 

T1 0 

T2 L°. 


(3.73) 


Equation (3.73) is expaned giving 


[XI -Ay] VyJ - [BFg] Vyg 55 0 


(3.74) 


[XI - E] v T - =* 0. 


(3.75) 


It is assumed that the eigenvalues for [Ay] and [E] are distinct 
since their locations are arbitrarily assigned as discussed in Chapter 


Suppose x is an eigenvalue of [Ay] but not [E]. Then [xi-E]”’ 
exists. Premulti plying equation (3.75) by [XI-E]’* yields the result 


V T2 * °- 


(3.76) 


Hence equation (3.74) is simplified to 




V 


[H-A t ] Vjj a 0. 


(3.77) 


Therefore Vyj Is an eigenvector of [Ay], 

Now suppose X Is an eigenvalue of [E] but not [Ay]. Then 
[AI-Ay] 1 exists. Equation (3.75) Implies that Vyg must be an eigen- 
vector of [E]. Premultiplying equation (3.74) by [AI-Ay] 1 and rear- 
ranging terms gives 


* n ‘ txI -\ r [s y V 


(3.78) 


Hence for X that are eigenvalues of [Ay], 




(3.79) 


where Vyj is an eigenvector of [Ay]. For X that are eigenvalues of 

cel 


[ (A I -At) B F 2 1 

. , J #T2 


(3.80) 


where Vyg is an eigenvector of [E]. Thus it is shown that a reduced- 
order model can be used to design a constant state feedback controller 
for a full-order system. The eigenvalues and eigenvectors assigned 
using the reduced-order model are retained in the full system while the 
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eigenvalues and eigenvectors not Included in the reduced-order model 
remain unchanged in the full-order system. It is also shown that a 
reduced-order observer can be used to Implement a full state feedback 
design without affecting the eigenvalues and eigenvectors of that 
design. 









CHAPTER 4 
OESIGN PROCEDURE 

. The methods described in Chapters 2 and 3 are the basis for devel- 
oping a design philosophy, and then a corresponding design procedure, 
for constant state feedback controller design. The design procedure • 
presented in this chapter is most useful when a designer is able to 
characterize the desired system in terms of the closed loop eigenvalues 
and eigenvectors as well as the time response. This chapter reviews an 
existing spectral assignment design philosophy. A corresponding design 
procedure and computer aided design package [11] are discussed next. 

Then an extension of the design philosophy is presented followed by a 
corresponding design procedure. This design procedure is included as a 
supplement to the computer aided design package. Lastly, the signifi- 
cant portions of the additional computer aided design software are de- 
scribed in detail. The modified design procedure uses reduced-order 
models and reduced-order observers with spectral assignment methods to 
reassign selected eigenvalues and eigenvectors in the full-order system 
model . 

4.1 Design Philosophy for Full-Order System Models 
The constant state feedback design philosophy for full-order system 
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models is illustrated in Figure 4.1. The objectives faced by a system 
designer are often many and sometimes conflicting in nature. However, 
the location of eigenvalues and eigenvectors, and the system time re- 
sponse are generally the prime consideration. After these objectives 
are satisfactorily achieved, secondary design objectives are considered. 
These secondary objectives include feedback gain reduction, minimization 
of closed-loop system sensitivity to modeling errors or parameter varia- 
tions, and noise suppression. The spectral assignment design procedure 
achieves a satisfactory control design by selecting an appropriate set 
of eigenvalues and approximating a desired set of corresponding eigen- 
vectors. The eigenvalues determine the rates of decay of the various 
system modes while the eigenvectors determine the relative contribution 
of each mode to the different system states and outputs. After a satis- 
factory time response is achieved with an initial eigenvalue and eigen- 
vector assignment, the secondary design objectives are considered. The 
freedom available to select the eigenvectors often provides a means to 
drastically improve these secondary design objectives while only slight- 
ly modifying the initial eigenvector assignment and thus the time re- 
sponse. This improvement is accomplished by modifying the eigenvectors 
within an area local to the original assignment. The direction and 
magnitude of the eigenvector modification is determined by a gradient 
search procedure as discussed in Section 2.6. 

4.2 Design Procedure for Full-Order System Models 
A computer aided design package written by Marefat [11] currently 
exists and is illustrated in Figure 4.2. The software package consists 
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of several special .purpose subroutines that are accessed by the main 
control program. The subroutines may be entered In any order to Imple- 
ment specific design objectives according to the design philosophy In 
Figure 4.1. The system description (l.e.. A, B, C) Is entered In Mode 
1. An arbitrary set of eigenvalues Is assigned In Mode 2, which then 
formulates the allowable subspaces for the eigenvectors. The desired 
eigenvectors are approximated in Mode 3 by projecting them Into the 
allowable eigenvector subspaces. Mode 4 allows the system designer to 
observe the time response for various Initial conditions and system 
inputs. The Initial design is then Improved by alternating between 
Modes 4, 5, 6, and 7 until a compromise between primary and secondary 
design objectives Is achieved. Modes 5, 6, and 7 modify eigenvector 
components, reduce feedback gain, and reduce system eigen-sensitivity, 
respectively, using gradient search procedures. 

4.3 Design Philosophy for Reduced-Order System Models and Observers 

The design procedure discussed In the preceding section is useful 
for systems where full state feedback is feasible. Another feature of 
the procedure is that It assigns all eigenvalue and eigenvector loca- 
tions. A system designer Is often satisfied with several open loop 
eigenvalue and eigenvector locations In a large system. The reassign- 
ment of the remaining eigenvalues and eigenvectors is better accomplish- 
ed using a reduced-order model that contains only those eigenvalues, due 
to reduced requirements for computer time and memory. Also, large syS- 
terns typically have fewer Independent outputs than states. A full state 



feedback. design 1$ Implemented In this case with an observer system. 

The observer estimates the system states In order to Implement the feed- 
back control law. Since some of the states can be obtained from the 
outputs, only the remaining states need to be estimated with an observ- 
er. A reduced-order observer Is desirable in this case. It Is designed 
using less computer resources than a full-order observer. Also, less 
hardware Is required for Implementation of a reduced-order observer. 

A design philosophy that uses reduced-order models and reduced- 
order observers Is Illustrated In Figure 4.3. In order to design a 
control system using spectral assignment with reduced-order models and 
reduced-order observers, a designer must have knowledge of a desired set 
of system eigenvalues and eigenvectors. The original open loop eigen- 
values and eigenvectors are compared with the desired eigenvalues and 
eigenvectors. A decision Is made as to which of the eigenvalues and 
associated eigenvectors are satisfactory and which need to be reassign- 
ed. The spectral assignment design procedure Is used to assign the. 
desired eigenvalues and approximate the desired partial eigenvector as- 
signment using the reduced-order model. Error between the Initial 
eigenvector assignment and the desired eigenvector assignment Is then 
reduced by a gradient search. Next, the resultant reduced-order feed- 
back matrix Is transformed to the full-order system. 

If all of the states are simultaneously available for measurement, 
then the full state feedback matrix Is Implemented. However, If some 
states are not available, then a reduced-order observer Is designed. 

The eigenvalues of the observer are assigned to be slightly more 



Figure 4.3. Reduced-Order Design Philosophy 
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negative than the dominant eigenvalues of the closed loop systen design. 
This Is done to ensure that the observer can respond quickly enough to 
follow the states being estimated. Theoretically If the observer eigen- 
values are assigned to be very large negative nunbers then the observer 
will provide a better estimate of the states. However, this Is not done 
In practice because the observer then acts like a differentiator and Is 
very susceptible to noise. 

This philosophy and the synthesis methods described In Chapter 3 
are used to develop an extended design procedure that exactly reassigns 
an arbitrary subset of the original system eigenvalues which are Includ- 
ed In a reduced-order model. A partial eigenvector assignment Is then 
approximated for these eigenvalues. This control Is Implemented with a 
reduced-order observer If there are fewer system outputs than states. A 
contribution of this thesis Is that the reduced-order design and imple- 
mentation are accomplished with the knowledge that the eigenvalues and 
eigenvectors not Included In the reduced-order model remain unchang- 
ed. 


4.4 Design Procedure for Reduced-Order System Models and Observers 
The computer aided design package written by Marefat has been modi- 
fied as Illustrated in Figure 4.4. An additional mode (Mode 9) has been 
added to Incorporate the use of reduced-order models and reduced-order 
observers In system design. The full-order system description Is enter- 
ed In Mode 1. If a reduced-order model Is to be used In the control . 
system design, Mode 9 is entered. Otherwise the design procedure con- 

















tlnues as described In Section 4.2. 

A flowchart for Mode 9 Is shown In Figure 4.5. Several of the 
existing program subroutines Including Modes 2, 3, 4, and 6 are called 
from Mode 9. Three new subroutines are also called from within Mode 9. 
These subroutines are described In Sections 4.5 and 4.6. 

After Mode 9 Is entered the reduced-order model Is formulated. 

Mode 2 Is then called automatically and the reduced-order model eigen- 
values are assigned. The designer Is* now prompted to enter the desired 
partial eigenvector assignment for the full-order system. An Initial 
eigenvector assignment is calculated from the reduced-order model using 
equation (3.24) and the assignment Is approximated using Mode 3. A 
gradient search Is then Initiated In order to decrease the error between 
the desired and actual partial eigenvector assignment. Upon completion 
of the gradient search, the reduced-order model feedback matrix Is cal- 
culated and transformed to the full-order system coordinates. A re- 
duced-order observer Is formulated next. Eigenvalues and eigenvectors 
are assigned to the observer using Modes 2 and 3. Finally a time re- 
sponse Is calculated and displayed for the combined system using Mode 4. 
If the designer Is not satisfied with the time response Mode 9 is re- 
entered. 

Two portions of the above design procedure required an extensive 
programming effort. Calculation of the cost function used In the gradi- 
ent search is described In Section 4.5 and the gradient matrix calcu- 
lation Is described in Section 4.6. 






Figure 4.5. Mode 9 ^ 

i 







4.5 Cost Function 


The cost function J that is used in the gradient search routine 
Is a measure of the error between the actual and desired partial eigen- 
vector assignments in the full system model. Calculation of the cost 
function Is accomplished In two parts. The actual partial eigenvector 
assignment is computed first using a subroutine called VACT. The actual 
partial eigenvectors are then used to compute the value of J in a sub- 
routine called ROCOST. 

A flowchart illustrating VACT Is given In Figure 4.6. The actual 
partial eigenvector assignment Is denoted by v^ and the eigenvector 
assigned in the reduced-order model is denoted by v^. This is consis- 
tent with the notation used in Chapter 3. The partial eigenvector as- 
signment v^ of the full order closed loop system that is obtained by 
assigning In the reduced-order model can only be determined after 
all reduced-order model eigenvalues and eigenvectors are assigned and 
the feedback matrix F is computed. The top half ®f equation (3.23) is 
given by 


v n- f v u ♦ v i2 cy-vViv (4.D 

For a real eigenvalue the subroutine computes v^ using equation 
(4.1). If the eigenvalue Is complex the calcuffiation becomes 
slightly more Involved. Separating equation (4.1) Unto real and imagin 
ary parts yields . • 
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(4.3) 


Equations (4.2) and (4.3) are used to compute partial eigenvector as- 
signments for complex eigenvalues. The actual eigenvector assignments 
are then used In subroutine ROCOST to compute the value of J. 

A flowchart illustrating ROCOST is given In Figure 4.7. If the 
desired partial eigenvector assignment is denoted by v Q and the actual 
partial eigenvector assignment is denoted by 7, then the cost function 
is calculated by 


J 


+ E 

U 


(Vij ’ VD ij )2 “ ij 


(4.4) 


where are arbitrary weighting constants. The weighting constants 

determine the relative penalty between the eigenvector component errors. 
For example, if one eigenvector component has a much larger weighting 
constant than the others, then an error in that component receives a 
much greater penalty than other component errors. 


4.6 Cost Function Gradient 

The cost function gradient matrix is computed in subroutine ROGRAD. 
A flowchart illustrating ROGRAD is given in Figure 4.8. It is seen from 
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equation (4.4) that the cost functon J Is a function of the partial 
eigenvector assignment 7. Hence, it is also a function of the desig- 
nator matrix X which is discussed in Section 2.6. By computing a 
gradient of the cost function J with respect to the elements of the 
designator matrix X^, it can be determined how to vary the designator 
matrix in order to reduce the cost function and therefore the error 
between v Q and v. Recalling equation (2.30), the gradient matrix is 
defined to be 
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To evaluate aVp^/aX^, q is substituted for i and the subscript 1 
is dropped In equation (4.1) to yield 


/ / I 


v q * t V ll* v i2^q I_A 2^ 


(4.8) 


th 


Since 7„„ is the p element of vector v then 3v„„/ax.< is the 
pq K q* pq ij 


p th element of 3v q/ 3 *jj* Solving for av^/a^j, then gives 


3v - ar v nVi i 3[Fv„] 

* 2l ^ Lsl * w*t> * “ 




ax 


U 


ax, 


ij 


a V gy ^ 

■ v u — * ,(» i-a 2 )'V][f - 2 -» — vj. 

jy AC q C C jsy q 

3X i j 3X ij 3X ij 


(4.9) 


In order to evaluate 3F/aX^, equation (2.73) is modified to be 


F ■ W 9" 1 


(4.10) 


f'h f u 

and the element of F in the p row and q column is denoted by 

f so that 
PQ 


f nn • R n [w][v“ l ] C . 
pq p'- q 


(4.11) 
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* row vector with a one in the p** 1 column and zero elsewhere and 
0^ is a colunn vector with a one in the q^ row and zero elsewhere, 
so that 


3f _ 

ax 


H 3 .R fjw_v-i + „iti] C 

v » u . .. J i 


8$-*- 


u 


ax 


1J 


ax ' <* 
3X ij 


(4.12) 
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Solving for aw/aXjj yields 


aw b __a_ 

8x u 


[«i. ...» 


(4.13) 
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From equation (2.68) , 


"j ■ -"j/j- 


(4.14) 
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Therefore 


aw 

ax 


c-m 

tj ‘ ax U x ,x lj 


(4.15) 


,th 


Noting that only the j column of X is dependent on X^, 


aw 

— * Co, .... o, r-M ] ax /ax o, .... 

3X ij 1 J J 


o]. 


(4.16) 


Furthermore, only the 1th row of Xj Is dependent on X^j. It Is 


,th 


easily seen that 3x j/ 3x jj a colunn vector with a one In the 1 row 
and zero elsewhere. Hence, equation (4.16) is written 
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ij 


[ 0 , ..., 0 , , 0 , . . . , 0 ] 
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(4.17) 


where denotes the 1 column of It is noted that 
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(4.19) 


Since only the j tfl eigenvector is a funtion of X^, it follows that 


a v 3 ^ i 

v * [0, ...» 0, — 1 , 0, .... 0]. 


ax 


ij 


ax 
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(4.20) 


Substituting from equation (2.77) gives 


a? 


ax 


ij 


* [0 0, [N x ) aXj/aXj., o 0] 


[o, ..., o, [n. o, ..., o], 

A j 1 


(4.21) 
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l u 

where [N. L denotes the 1 column of [N, ]. To evaluate 3v„/3X<<, 

Aj i a ^ q ij 

equation (4.21) Is postmul tipi led by C^. Similarly, to evaluate 
3 v /3 X. , equation (4.9) Is premultiplied by R . 

KM K 

Hence, It is shown that the partial derivatives are computed by 
selecting appropriate rows and columns from the [N^] and [M. ] matri- 
ces. Tills reduces the calculation of the gradient matrix [GR] to a 
bookkeeping operation easily Implemented in a computer program. It Is 
not necessary to numerically approximate a derivative quantity. Sub- 
routine ROGRAD computes the cost function gradient matrix using this 
procedure. 
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CHAPTER 5 
DESIGN EXAMPLE 

The design procedure described in Chapter 4 is Illustrated in this 
chapter by an actual design problem. A controller Is designed for the 
lateral axis model of an L-1011 aircraft using a reduced-order model and 
a reduced-order observer. The resulting design is then compared to an 
output feedback controller designed by Andrey et al. [16]. It is shown 
that the design procedure presented in this thesis is a viable tool for 
constant feedback controller design, 

5.1 Original Lateral Axis Model 

The lateral axis model of an L-1011 aircraft Is used as the orig- 
inal full-order system model. The state vector x is qiven by: 

xi * r * Yaw rate (Radi ans/second) 
x 2 * 6 3 Sideslip angle (radians) 

X 3 » p * Roll rate (radians/second) 
x 4 » $ 3 Bank angle (radians) 

xs 3 $ 3 Rudder deflection (radians) 

x 6 3 3 Aileron deflection (radians) 

x 7 = f = Washout filter state. 

' w 
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Rudder and aileron deflections (states 5 and 5) produce changes in the 
yaw rate, sideslip angle, roll rate, and bank angle (states 1-4). The 
coordinate system is Illustrated in Figure 5.1. Under certain con- 
ditions yaw rate is equal to the derivative of the sideslip angle with 
respect to time while roll rate is equal to the derivative of the bank 
angle with respect to time. The washout filter is a high pass filter 
for the yaw rate. 

The A, B, and C system matrices are given by 


*0.154 

*0.996 

0.249 

0 

0 

0 

0.5 


1.54 

-0.117 

-5.2 

0 

0 

0 

0 


-0.0042 

-0.000295 

- 1.0 

1.0 

0 

0 

0 


0 

0.0386 

0 

0 

0 

0 

0 


-0.744 

0.02 

0.337 

0 

- 20.0 

0 

0 


-0.032 

0 

- 1.12 

0 

* 0 
-25.0 
0 


0 0 0 0 0 
0 10 0 0 
1 0 0 0 0 
0 0 10 0 


The system Input, u, consists of components 




80 


u i * s rc * R*‘ ;, der command (radians) 


u 2 * 6 ■ Aileron command (radians) 


The open loop eigenvalues of this system are: 


A 1,2 * -0.08819 ± j 1.269 - Dutch roll mode 


-1.085 

-0.00965 

- 20.0 

-25.0 


- Roll subsistence mode 

- Spiral mode 

- Rudder mode 

- Aileron mode 

- Washout filter mode. 


The open loop system time response Is shown In Figures 5. 2-5. 8 for zero 
Input and an Initial condition of 4>(0) = 1 degree. After ten seconds 
the system states are still oscillating and the bank angle <f> has not 
yet reached zero degrees. 

It Is known that a desirable eigenvalue assignment for the system 


X l>2 * -1.5 ± j 1.5 


A 3 , 4 * -2.0 ± j 1.0. 


When the roll subsistence mode x 3 and the spiral mode x 4 are a com-. 
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Figure 5.3. Sideslip Angle-Open Loop Response for <j>(0) = 1* 
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Figure 5.4. Roll Rate-Open Loop Response for <)>(0) = 1* 
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plex conjugate pair they are collectively referred to as the roll mode. 
It Is also known to be desirable for the roll and dutch roll modes to be 
decoupled. This decoupling Is accomplished by the eigenvector selec- 
tion: 





where X denotes "don't care."* /Indry, Shapiro and Chung [15] closely 
approximate the above eigenvalue and eigenvector assignment for this 
system using constant output feedback. Eigenvalue/eigenvector assign- 
ment techniques are used to design the constant output feedback matrix 

j 3.35 -0.159 -4.88 -0.379 

* * 1.42 2.38 -6.36 3.8 

The closed loop time response is shown in Figures 5.9-5.15. The closed 
loop eigenvalues of the design are 

* 1,2 “ -1.052 * j 1.497 

X 3 , 4 « -2.001 ± j 0.9995 

* -17.05 
X 6 - -22.01 

X 7 * -0.6989. 
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Figure 5.12. Bank Angle-First Closed Loop Response for $(0) * 1* 
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Figure 5.13. Rudder Deflection-First Closed Loop Response for $(0) = 1* 
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Figure 5.14. Aileron Deflection-First ri * co m i 
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Figure 5.15. Washout Filter-First Closed Loop Response for $(0) * 1* 
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The first four components of the first four eigenvectors ‘are: 


' 1,2 


0.03066 ± j 0.3488 
-0.0036 ± j 0.0004 
0.0013 ± j 0.0011 


'3,4 


-0.0029 + j 0.0012 
0.0045 ± j 0.0053 
1 

-0.3999 + j 0.2000 


It is noted that by using constant output feedback that the four eigen- 
values X x - x 4 are placed almost exactly and that the roll and dutch 

roll modes are decoupled. However, the other eigenvalues (As - X 7 ) are 
also moved by the design. The design procedure described in Chapter 4 
is now used to formulate an alternate design that exactly places Xi - x 4 
without changing x 5 - x 7 . Eigenvectors for - x 4 are also assigned 
to achieve roll and dutch roll mode decoupling. This is achieved with- 
out modifying the eigenvectors associated with x 5 - x 7 . Furthermore, 
this design is done using a reduced-order model to specify a constant 
feedback matrix for the original full-order system. The full state 
feedback matrix is then implemented by dynamic output feedback. 

5.2 Reduced-Order Model Design 

Since the open loop values of x 5 - x 7 are known to be desirable, 
and the rudder, aileron, and washout filter states are unspecified, no 
reassignment of these modes will be made. Therefore, they are not in- 
cluded in the reduced-order model. On the other hand, X x - x 4 are to 
be reassigned and are included in the reduced-order model. The full 
order system matrices are transformed by equation (3.3) and partitioned 


A 
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as In equation (3.5) to yield the reduced-ordeh model system matrices 


* 



• 

• 

m 

■ - 

-0.08819 

1.269 

0 

0 

* 

-1.706 

-0.02580 

-1.269 

-0.08819 

0 

0 

. Bl* 

0.3961 

0.06988 

0 

0 

-1.085 

0 

-0.2772 

-0.2878 

0 

0 

0 

-0.009165 


-0.2698 

-0.1528 


Spectral assignment synthesis methods are then* used to assign the eigen 
v al ues 


* 1,2 * -1.5 ± j 1.5 
and 


*3 >4 a "2.0 ± j 1.0. 


A partial eigenvector assignment that achieves roll and dutch roll mode 
decoupling is given by 


" - 



20 


0 

6 ± j 7 


0 

0 

* V 3,4 “ 

20 

0 

-8 + j 4 

• 


. . 


An initial attempt is made to assign the eigenvectors using equation 
(3.24). The partial assignment in the full order system is found to be 


• « 


- 

19.44 ± j 0.33 


-0.10 + j 0.04 

6.76 ± j 7.25 

If _ 

0.08 ± j 0.10 

1.06 ? j 0.10 

* V 3,4 

20.20 ± j 0.11 

-0.42 ± j 0.64 

-8.54 ? j 3.91 
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The gradient search routine described in Section 2.6 is now used to 
improve the initial vector assignment. Elements of a weighting matrix 
are entered into the computer and a value is calculated for the cost 
fucntion J using equation (4.4). A cost function gradient is calcu- 
lated as in Section 4.6 and the initial eigenvector assignment is varied 
to reduce the cost function. The weighting matrix is varied to increase 
or decrease the relative Importance of each eigenvector component and 
the gradient search is continued. This procedure is repeated until a 
satisfactory improvement of the initial assignment is achieved. In this 
example the final partial eigenvector assignment in the full system 
model is given by 


- 


" * 

19.45 ± j 0.34 


-0.10 + j 0.04 

6.76 t j 7.25 

W - 

0.08 ± j 0.10 

0.45 ± j 0.33 

’ V 3,4 

20.20 ± j 0.11 

-0.07 ± j 0.68 

-8.54 7 j 3.91 


The vectors are scaled to give 


1 


-0.005 7 j 0.002 

0.35 ± j 0.37 
0.02 ± j 0.02 
0.003 ± j 0.04 

’ V 3,4 = 

0.004 ± j 0.005 
1 

-0.424 7 j 0.191 


It is seen from the above vectors that the roll and dutch roll modes 
have been significantly decoupled. The required gain matrix in the 
reduced-order model is given by 



\ 


C . r 1.319 -1.650 0.169 -1.7241 

[-3.854 °- 583 *5.810 31.18 J* 

The constant state feedback matrix In the full order system Is computed 
using equations (3.9) and (3.10) to be 


[3.66 

-3.13 

-0.176 

-0.372 

-0.137 

0.003 Ol 

|L54 

-6.03 

2.70 

4.35 

-0.008 

-0.120 Oj 


In order to Implement this full state feedback matrix, an observer Is 
now designed by the procedure described In Section 3.4. The observer 
eigenvalues, are selected so that 

*01 “ " 5 » *02 “ “ 6 » *03 “ “ 7 * 

This selection makes the observer modes faster than the modes contained 
In the reduced-order model. The observer eigenvectors, Vq^, are arbi- 
trarily assigned to be 



p» <m 
1 


m m 
0 


0 

V 01 ■ 

0 

0 

• V 02- 

1 

0 

U 

> V 03‘ 

0 

1 


The observer matrices are then calculated to be: 



100 


and 


It 


1.287 

-0.029 

-28.3 

-0.372 

1.184 

0.710 

-17.1 

4.35 

-89.2 

2.35 

-17.8 

0.124 

-42.2 

-83.7 

57.6 

-0.116 

-5.29 

• 

-0.035 

48.5 

0.270 

tB • 

[20 

°1 


- o 

25 . 




The observer Is now used to Implement the full order system feedback 
matrix F. The closed loop time response is shown in Figures 5.16-5.25. 
It is shown that the response of the yaw rate and sideslip angle for 
this design are more desirable than for the previous design since there 
Is less disturbance and faster settling time for both states. On the 
other hand, the roll rate and bank angle responses are almost Identical 
for both designs. The controlling surfaces and washout filter states 
are all well within the physical limitations of the system. This illus- 
trates that a viable constant full state feedback control system can be 
designed for a large system using a reduced-order model and that this 
feedback design can be successfully implemented with a reduced-order 
observer. 

5.3 Summary and Concluding Remarks 
A new design procedure for the control of large systems using re- 
duced-order models, reduced-order observers, and spectral assignment 
techniques is presented. A reduced-order model is formulated containing 
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eigenvalues of Interest from an original full-order system. A constant 
state feedback matrix Is designed for the reduced-order model that, when 
Implemented about the full-order system, reassigns the eigenvalues con- 
tained in the reduced-order model while those eigenvalues not included 
in the reduced-order model are retained in the full-order system. It is 
then shown that the full state constant feedback matrix for the original 
full-order system is implemented by a reduced-order observer if all of 
the system states are not simultaneously available. 
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APPENOICES 


Appendix A contains a software listing of the modifi- 
cations to the spectral assignment computer aided design 
package discussed In Chapter 4. This Is followed by an 
example of an Interactive design session In Appendix B. 
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i:ml wux rr< vo ihv. a. ms. ns. ns.io.io.v. lu.irfei 

Wl vnm.f'l IU.VO.HS.HS.NS.IO.IO.AI AM. 10. ICR) 

Ml Al« I* I. MSS 

DU «I4 J-l .NS 

APSAl .ADS1A1 AMII.JII 

«l < ADSI 4. .L r . ZEXOIALAHI I , JI-fLUAI 10) 

ai a cuwmhof 


c paosl s 

CML USUrmSWIALAM! •S.ALAM, 10.MS.MT..4) 

CAII lmULM’IVniMV>«,MS,KS,KIclD,10,M^M,IDtIER> 

c pmiu: 3 

C CALL USI*'M(SISMJIMI>S.M.MI l ia,NS,MI,4> 

CAU VIKSLI r(C,Vn,M),MS>KS,IO>tO,aAM,lP,tElt> 

c rAitsj a 

C CALL US.M'H(5««CLAHI,3>Ct AM. 10. MO. MS. 41 

CMfltltltmt PAR) HUM TILDA SYSTEM 
Ml AOO 1*1. MO 
Ml A60 

MAMlHI*J)>MiWI(t»J) 


AOO CONTIMK 

Ml SVV 1-1. MS 

do syy j-i ,ms 

IK <1.07. NO I .AND. « J.BT.MO) IAUK27(I-«ViJHnt>fLAMa,J) 

s.y cttmmsc 

DO A01 I - 1 . KO 
Ml AOI J-l.MI 
DLAMKI.JI-SLAMCItJ) 

Aei ruMimiE 

Ml AOS I- I .MO 
Ml AOS J-fttMO 
CLAMKI.JI-CLAMd.JI 
AOS CtUlTIME 

c**»»«e»«9*08* DISPLAY RO MODEL ttMlIttmUtMUMOMMItUNMIM 

URITE«A>004) 

804 rnfcnATliX.SSi&IISI! TO BISPLAY HO MODEL Tl 
Kf ADIS. .IKK 
iriKK.LE.OIOO TO BOS 
SMI Tl'(«»A03) 

403 FMtnAI(IX.///>SOH REDUCED ORDER MODEL-//) 

CALL IISM'HCVSMATRIX Al.y.AIAMH.lO.Mfl.MO.41 
CALL IROrMIVIMATRIX »1.V.D4 AMI. 10. MO. MI. 41 
(SAIL IRSSPMI VI MATRIX Cl.y.CLAMl.lO.MO.MI.A) 

ciiimmm rt.KPORARiLr store orictnal stbiem immnhmmimmii 

DOS SR I TEI. lS.REC-SIMS.nl • MO 

MRtTE(.IS.REG>3>«<Aai.IJ).IJ-I.MSI.It-I.ND) 

HRm:<32.RL'i:-4M(B<IItIJ>tIJ-I.Ml>.II-l.HS) 

URDEIIS.REC-SI IICII I. I J»rIJ»l.l«i>.lI-l.llO) 

»•*(***** STORE RO M01A1. INTO FILE *SO 



r~ . 

in 
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# 


in* \w\t\-: 


.K 


' \ 






f r 

v~' 

i 


mr Tr<^p.irrc-i iMa.ffi # md« f mi «znto 

UN I T E I i» 0 . Rf C * 2 > I I At AMI I I I 1 . 1 J > , I J* ! • MO > • ! I - 1 • NO > 

hr ! n:< ;*»*Ktrt-;i x i m.mu « 1 1 • i j> • i j*t »mi > • 1 1-1 • no> 

UNI 1EI20.RCC-4 I < (CL AMI C I I • I J> « IJ-I .MO) • I 1*1 .NO) 

Lttfttl IK) TO HUPf 2 Ml* <PSIiiN CIWKVAUAT- FUR RO SYSTEM $$9999$$ 

UK 1 n<6»A04> 

AOA TORIMI <lX.iRMH*>.3AN REDUCED UNDER EIUMIVAt 1C ASSlOMttf NT * t /(IH4 > 

i.vtu Hftnt:? 

C$999999 INTER Dtf.IMD t-'ANT 1AI CIUf-NVCClOR AV41(U»2NI MllHMl 
UNI1CIA.AOS) 

A<*3 I (IRMA Til X.44ICN TER DESIRED EAR MAE Ellir.HVLClIIR AGUlURCNfl 
.1-4 

AOV CXMEIA.AOAhl 

AOA KWHA;! JX»UH4 iw.wvr.rnm V»I2> 

readis- • > i vici. < i • .f > • vr< i > , i -t rW) 

it <ARH<t iN4.i)i.Lr«/rRii)m) to ado 

DO 60 / 1*1 »Ntl 
VICSI I. -Ill )-VHI> 

AO/ lONTlROE 
XIH 
AOt« J*JU 

ir cj.le.hu >imi ro aov 

C999999999 C/V tU Ml. AMO DISTLAY INITIAL IRJTS8 Itmilllltltllttllll 
mi A10 l-4.NO 
IK* AIO J"t.NO 

aio nvMfiMur 

CALL LIKV2F CU.WOrlO.AHAT.ID8T.yKANEA.ItR> 

C CAUSE 5 

C CALL ICiW'NI/HVUIMVI./.IUIAT.lO.HO.MO.A) 

CM L VWU FTIAHAT. VBC8.N0.NU.N0. 10.10.8. 10*100 
URXTEIA.AIl) 

dit from* rux •*?}«** t« riunnm v matrix for initial, msimcrti/i 
iv 9 h wrt*:t%*£* which v arc complex) 

CHARI 1 >** INI TtAL GO* 

CHARI?)"* CSS FOR VI* 

CALL IJSW*«<CHAR.2O.G.IO.«0.MO.4> 

Clttllttltm ASH I UN INITIAL EIGENVECTOR BUCNI •IMttUlllllltmitl 

FAI L modes 
IV* 20 

NT AO IIU.REC-l) N&.N1.NQ.IK3T.ZCRO 

WAD IIU.REC-2) IIAIII.IJ). IJ*1 .MS) » 11*1 *N8) 

READ <IUtRCC*3) I.H4rNIKII-l.Na> 

RCA0I JU*REt~4> ( (Cf II* IJ)*I>I.m)iIM .NO) 

C rAUSC A 

c uRjrriA.eio) 

•JO FORMAT I IX.3/HTHT FQLLOUIN0 18 RO MODEL AFTER MODE!) 

C CAiL IISUFNI2HA1 .2.A.tO.WS.N3*4) 

C CAIL IIS4N f1(2ft9l .2.8. 10.N8.N7.4) 

C TAiiSC '/ 

C CALL USWFH(2MCI.2.C.10.ND.|tS*4> 

HO JO .Hit NS 
IU-20IJ 

ir#.N c access* '* iRrcr'»RrcL"?02 

l.DNIT-UO 

RtAD <IIJ»RrC«l) LRCCJI.LIHU) 

10 CUNflMK 

C ' PAUSt H 

. Wl !M2 1*1. MS 

C WfniKIA.RlDI »LRE! I ) .LIMII) 

011 FORMAT! IX. All.AHDBA. 12. 3X.3HREAL*iFA.3f3X.3HIHAn*.F4*3) 


HI? CUNTIWJE 

IUT*20JNStl 

trCN I F II. E-'CtSOtNI * .ACCESS* 'DIRECT ' *RL(L"102 
i.iaMir*iur> 

read uur.Rroi) C(Viu.tj).ij-i.iif;>.ii*i»Mr.) 

REAP !IUT.ktC*2) (CXXdl.IJI.IJ'f .RS).II*J.KI) 

WAD (UIT.RE04) CfFf II.IJ)* IJ*1.M&). II-l .NT ) 

RtAP (IIITrM'C^) KAMA? (II. IJ)»IJ^t.NS>. 1 1*4. NT) 
f. FAII IISUEM MOWV.fklX V I. 10. V. 10. NT.. MS. t) 

L CALI (KWril I lOHrWifKIX XXI . 10. XX. 10.NI .NS *4) 

C PAUSE V 

CHARI | >**UAIM MATkl* 

MIARI2)«*K FI* 

CAI.L USUFM (CHAR* 14. F « 10. Ml »NS.4> 

C PAUSE lO 

MIAMI I >**MATRIX aha* 

MIAMI?)** T I * 

C CAIL tl&UEM 4 CHAR. 12.AHAT* 10.NS.NS.4) 

REAP<J?,kEI>2>URIUO.NlX.NOX « 

X2-0MIMI NS 

C$99999999 ASSIGN ICIUHTIM! CONSTANTS FUN OCAOILIfT SEARCH Ittlltl 
90S DO /O I*i. NS 
DO 70 .HI. NS 
ALII. J)*F1 OATH) 

TO COMTINVE 

URMCCA.Vy) 

99 FORMA TC1X.20C IMS) .2/H RO CICMCVCrTOC inrR0VCHFNT.2OIlH4)) 

UR)TC(A.F) 

9 FORMAT! IX. 20H UK) (SIT I WO CONSTANTS) 

CALL VSUF(1imiUCieKrSI»I.AL.I0*N8.N8.4) 

UUllLIA.lt) 

U FORMAT! IX. 1AM UJSH TO CHMBEf) 

RCAD IS.t) KK 
IF IKK.LC.O) 00 TO 30 
UNITE (A. 3) 

3 rDKNAf ! IX. 17HTNfCR NED VALUES 1 ) 

REAR IS. 4) !!AL(I.J).^*|*NS).I«l.|f9) 

ciMiiiM cwswri oraricnt SEARCH itietmtimiminttititmit 

30 CAL L ROOCSTICJ) 

URITL !A*4) CJ 
4 FORMAT !f X«SHC0ST«.C19.4> 

FAIL ROOcAD 

CAIL HLARCHlCJ.ROCttSI.ROCItAD.S) 

URITL (IUT.REC*!) (!V!II.IJ).I>*t.N 3 >» 22 ««.NS) 

URITC (IVT.REC-2) € T ir nJT . 1 7= 1 ; gg J .1 1- r, ET) 

UNITE !ItlT.RrC*4) I (Fill. I J) t I J*1.R3).II*1.NI) 

UNITE !IUT.REC*3) (<AMAT!II« KJ).I>1 »NS) .11*1 .NS) 

cakl usurn ciomhatrix v i.io.v.io.ns.ns.a) 

CALL VALT 

CAIL UHFKUMUrfttX WS.IO.VFS. I0.NS.NS.4) 

URITEI A.R04) 

•04 FC9VtAT!IX.24MUANT TO CONTINUE SCARCttt) 

RFAOIS. $)KK 
II IkX.t)T,0)GO TO TOS 
C MMITI IA.902) 

902 FORMAT T1X.44MVISH TO DISTLAT THE NORMALISED ETCCNVECTORST) 

L REAP (5.4) KS 

C IF IK8.LE.0) 00 TO 903 

o> 


t 




Vi 


a 


4k.- 
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C CALL DSCLAY LNS./CRO) 

*03 CONTINUE 

C l ALL IISMTH <10flMAr*lX XX|'lP.XX»10»lt!»NS»4> 

CllARL 1>»* GAIN MA* 

nWK«.*!l-*TMX ri* 

cail uswfh cD<Aft*i4»r*io*iu*M$»A) 

f'HAttt 1 AtIA* 

CALL Fll^rtAJtAI'NS.lO.Z.UEia.Z.lO.ircAREAtlER) 

c call ti!rav<w&AM~WA9«a*iff:ta'its«i«4> 

CUAk<2>-HI* 

call usurn <cuAR.i:»*AiiAr*io.tfs*Ns.4) . . 

cm*«m accemd zeros ro feedback matrix ana xr«m 10 original cot 

»0 701 I»i *NI 
DO /OI J'NOU.URIGO 
rci#j>*Fi.OAfio> 

REAM J2.RI tM)( <01 IX* I J) • IJ"»1 tNllt II-ltKSI 
REAmJ2.fttT«5)<CC<fl9lJ)f i>-t *ksi*ii-i *m»i 
Kt All I 32. RF 04 > << r m • IJ) 9 1 J- 1 .NS I v 1 1 •! till > 

cail iJsurrt<5Hr-Arr**s»F«tO'iti'f?s»4> 
ctimni cunrmx wuwi matrix n i hikv miitmitmtii 
DO 703 1*1 *NS 
DO 703 J>»i*NS 
M< I# J)-*FLPAT<0> 
ira.ro.j)H(i»j>«FLOAm> 

7rO.LE.NH)K(l>JX(lfJ} 

703 CWUIWX 

CALL LIKV2F (H.US* lOtMINV* IPOI.UKAREA.IER) 

C CAUSE 13 

C CALL USUFn(SHrtlNVIt3*mUV«tO»nSfNS*4> 

C CAL L IfSUFHT 21*11 *2*H*10*N3*NS*4) 

ctmtttt xrnRfi system ako f matrix nmtmmimttmi 

CALL VKtflLFf- <M'A,tt3fN3iNS»10«t6'AfrMT»10»U.*> 

CALL V?3Jt.Fr<AHAf *HIUV*ftS«MS*ISS* 10# 10.ALAI1* 10* ICR) 

CALL VWXrFLM.B,KS*US*WI*I0*l©»*CAM»10*lC») 

CALL VfflJLrr TC.MINV.IM.MS.HS. 10# 10. CLAM# 10# IOU 
CAt.L VMtJLrF(F#MIir7#NI*NS*NS*f0#l0#AKAT#10#IERI 
no his i*i«ki 

DO 840 >liNS 
»LAmifJ)»AHAT<I*J> 

840 ctmitfie 

cmtltttt DISI1. AY XrOftMFO SYSTEM Itmtllttlllttl 

WRI1ET4.H01) 

801 FORMAT < 1X»30H UISH TO DISPLAY TILDA SYSTEM?) 

ki:aim8»dkk 

If tkK.LE.OHKI TO 802 
URITC<4«704) 

704 FORMAT! IX* 2011 SYSTEM XFORMED »T «*//> 

C I’AIISt 15 

CALL ffSUFff < /HA rUDA*7*AlAM*10*ll$«MS#4) 

CALL !l?;yrM<7MP TILDA. 7*Bt AM* 10.MS.fl) #1) 

CAUSE 14 

call tmurm/HC rii.DA«7#a.AM#io»im»ies«4> 

CAIL UflLTMCTHF Til DA* 7* AtfAT * 10.ftI.NS* 4) 

CmtttlttMttt CAR Till {XI TIC SYSTEM tmilttttti 

802 Wl 70% 

Ml 70S JM.WJ 
ALAMI Jl !».f>~ALAM<f #J) 

ir< n tin. i;r.L*s>mi ro 823 

Al AM?1<I*.I> 'At.AMW IMIltJ) 


701 CONUHtJF 

C PAUSE n 

C CALL IISUrN(10HrArrLUlCD!*10.FtlO.Nl*URI0O*4) 

tall mirr Tr*voir/v»MT*aRrrm*aRiix)#io.to*Aiu.T#io.zER> 

DO 702 1*1*111 
00 702 JM.nHIWI 
rci .JWJIAK J..I) 

702 cor; nine 

CAUSE 12 

CALL IJSiriUIOMr XFURMEDt #1P»F« 10*NI»0RIGT)*4) 

UNITt <32>RrOAH<r<Il9lJ)9l>*l *(*1001.11*1 .Mil 
C CALL (*UrM<2ICI«2#r#10#l<l*<ttI8O#4> 

Cltllltll RETRIEVE ORIGINAL SYSTEM DATA ft STORE F 
4 TV RTAD(32*REC*2)KS«lir*IIO 

READ( 32.RTC *3) < LAC II » I J)« I J*1 «RS> t tl*l #KS) 

823 irCJINO.|IT.ff*;)t«) TO 824 
A9AMI2<I#J>«AI AMTI.JfNO) 

024 IT CLIIMn.OT.NSKAttD.f JIMn«0T*!tS)>8O TO 705 


AC.AH22C I » J)»A| AM ('I HOt JUKI) 

705 CONTINUE 
C CAUSE 40 

C CAI L USUFMC 7HAT1L 121 • 7.AL AMI2.10.MO* 12*4) 

C CALL USStS'MC 7HA f IL22I # 7* ALAM22* 10* 12* I2«4) 

DO 704 1*1. NO 

DO 704 J™1 *KI 

BLAM1 <1 #J)*BL AML 1#JI 

IF < II m .LE . MS ) DLAM2 ( I « J> -8LAMC I IRn. J) 

704 CCHTINME 

DO 707 I* 1 * MU 

DO 707 J^i.Nft — 

CL AMI C I * J >-CLAM< I * Jl 

If < J I NO . LE . ie; IU.AM2C 1 9 J) -CLAM ( f » JIMtl) 

707 town IRC 

DU 70S 1*1. MI 
DU 70S J~1.N0 
I TILUI* JI^AUATCI #J) 

I F C Jf NO . LC • NS IF T 1L2 L 1 9 JI-AHATCt* JUfO) 

CONTINUE 

* ’** W ” 2 T0 ^ i8I0 “ OWWWOt BYKMUCS 

•» rOf-«Ar<IK‘.MMS»,I2.3x,3HKO-,I2» 

MCDG-KOMO 

bo /ov i-i.m 
do /©» 

M AflI2IIrJ 
TO? CIMTINIJF 
C CAUSE l« 

C CAM. IKW*M«W:!-A»Afll2l.«,ALAMI?.IO.MO.W>»3,4» 

C CAIL USUrni7liALA«;2l.7.AI.AH22. 10, 110*3. MODS. 4> 

CAI.L IRAMS1 (MAH22.KnkS.Ktl9S) 

CAM TRAMSt(ALAflJ2.MO,MnftS> 

C TAIKE 17 

C CALL IK4»~HIVH-ALMU2TI,r.ALAM2flO,MttS.HD,4> 

C CAIL «»S'N(H<ALAM22TI.a.ALAM22.1P.Ma><S.M0DS.4> 

IFIKn.LE.WIBSMRI TO B31 

inr-i 


f 

I 

*• 

‘ j 



■g O 

09 , 
c r» 
x> r*. 


<C 

c *; 

C 

-i«a 







C*u .?»ALAM2» 10.KO8S.N0.4) 

8RITE<4.833)Nt>8S 

B34 r(«nAmX'21HY0UM*T 6C1KT MH1CH.I2.42H OUTPUTS IIIU 8£ USE* 

iruo tin obs^stoh 

IKITi <4.8341 

834 FORMAT! 1X.43HTHE OUTPUTS CORRESm® TO COLUMNS IK -At Mil 21) 

837 NUT TC! 4. 8351 

833 rCkMAf<IX,«OHSn.CCIED COLUMNS WILL CORN A RA1RIX THAI MUST 8E K 
1XNOUI AR.//.34HE»m*R 0UTPU18 TO 8E USED ( INTEUCRI) 

REAM 3 . * > II! R< K I » 1*1 * HUM) 

Wl 834 .f*l.NOBS 
DO 834 I»l»Ntl83 
AllAft I,J>-ALAM2(I.ILR< JH 

834 CUII f INI IE 

CAIL tCIWI rt(4IW nODl.A.AHAT.lO.NORS.KOftS.A) 

ICfc»0 

CAU. LINV2F < AJIAT .NOBS. 10* At • IDGT.ULAWrA. ICR) 

IF tlLR.tO* 12YHMJ TO 837 
WJX-MOPS 
00 TO IMS 
831 NUX-NO 

l>0 838 I -1. NUBS 
DO 1138 J-l.NO 
MIAra>J)»MMU3UfJ) 

838 CUUIJNUE 
843 CONTINUE 

Ctmm STORE ALAN22 t At ANl 2 A3 A ft 8 AND CHANGE Nfl Mtlltllt 
WRITE! 20. Rt*t>l INOIft.HIlX.NOt'S# IPGT.ZERO 
MR X TE ( 20 . rax-2 1 << AL AH22t 1 1 . Ml . I J « i .MOW; ), 1 1 • t ,w*D3> 

WRUF ‘ 20. REt>3 ) C C ANAT< 1 1 » 13) • W* t »WIX ) . 1 1* i .NOBS) 

W««( A;»!» I ON KIUf.WAI.tHTS AND C I UIHVEL lOK'i FOR OWaRVCR Mtltllt 

a t»st < tain - iut . status* *»xa w * 

00 815 1-i.HO 
.1*1120 

CLOSE !IJI1IT*J> 

8i5 continue 

UR I 1E!4« 710) 

71 0 r0RHAmX.20!lllt>.32H ASSIdtt CIOTKVALIVS FOR OBSERVER* 10( IHt 
CAI L nude? 

WkI1£(A./ll> 

/II FORMAT C IX. 13C1H8I.3JM ASSIGN ElOT.NVTCfORS FOR 8KSTRVFR* MUM 
CAIL nniUTJ 
WRITE (A. 7.13) 

7js roRNAT<iXr33ininoi ro reduce bain for odsirvcryi 
m:ai>(3»b>kk 
IF(KK.LE.O>QO III /34 

call mnE4 

/34 CONTINUE 

Cllimi RF. TRANSPOSE AIAN2? 8 -ALAN 12 t XPOSE F TO INI L Itttmitl 
MIX-ORZOU-KOBS 
IKI /12 1*1 # MODS 
W» /!2 J**l »NOX 
ALAMJ2<I».J>— A1AMI2U.JI 
712 CUNT INI IK 

CAl.l. TfcAMSl(ALAM22.NS.N8> 

CALI. TRANSMALAI112.N088.N0X) 

C PAUSE l» 

C FAI L USMFm7HALAH17| • 7. At AMI7. 10.K0X.MO8S.4) 

C CAI.L USUFNC7HAI.AN22I . 7.ALAH22.10.NQBS.H083.4) 


IUT*..\MNfiU 

KiMKiu7tn:cMmr(iZfij)ii>i»i»)rii«i*Nn 

irciFtt •ME.neo 10 «mi 

DO 83? I-l.HOX 
DO 83? >l»NOM 
MI *J) '"FLOAT <01 
83? cowmoc 

DO 040 I-l.NS 
DO 84C 3-1. NS 
L(ItR(n>JK(!»J> 

840 CONTINUE 
(XI TO 844 

841 Ml 842 1*1 .NUX 
Wl 847 J~1.N088 
tu»j)-F<r»J) 

842 CGHflKlS 
844 CONTINUE 

CLOSE ((Wl T*IUT. STATUS* 'KEEP' ) 

CALL TRAMS 1 I L. NUX. NS > 

PAUSE 1? 

XFLT -O 

LAI I USUFMC2ML I »2*l.« 10.WS.NDX.4I 

ctiittumut court i tc n MimiiiwiimimitliWM 

lit* At* C 32 . R£t> 2 ) MS . «1 • «0 
NOS S*M?5-NU 

CALL VtHtrMFUtl.L.KI .MODS* NO. 10. lfr'AMU' i»t IER) 
Ml 713 MI.NI 
Ml 713 .1*1 .NO 

C R IS RTPRCSEMTC8 8T FTltl 

FI IL1 1 I#.l)"FTItl<I.J)4AHAfCI# J> 

713 CONI IMN 


cmittmtim comiTE r tmttmmtm 


CAIl VfltATf (I . AI.ANI 2.N0BS. NO. NODS. 10. 10. AHA r»lO*KERl 
DO 714 1-1. WOOS 
DO 714 J^I.NODS 

CODS f I . J> -AI.AH22? I . Jl -ANA) < I • 3) 

714 CONTINUE 

ctuimtm courtm o tutMtmuimtmtii 

CAI.L VMULFF <FO*S.L.tflI*S.MODS.NO.IO« 10. AHA 1 . 10. UR) 
CAU VttXJ r r<l .At AH 1 1 .MU8S.N0.M1. 10. 10.AL AH27. 10. IER) 
DO 7*3 1*1. NODS 
DO 783 >> l.Nfi 

C O IS RCPRITSf.N TtD DY A1-AH2? 

M AH22 < I • JI*A1 AH21 C I . Jl -ALAH27U . Jl 4AHAT < I . J> 

715 CUNTIMUC 

CIMlimm C0HP1ITE 8- TILDA tltltlHIM 

' CAU. VNULMCL.BLAftl .NODS.KO.Nl.lO* 10.AHAT.10. IFR) 

DO 714 1*8. NOBS 
HO 714 JM.NI 

C AHAT Rtl'RFSENTS T-8rit.HA 

AHA T < 1 • .1* -8LAM2I I * J>-AHAT( I . Jl 

714 continue 

Ctllttimi DISPLAY R.8.E.8 18 

CALL USUFIKYHnATRTX El .Y.E0D3. l0.«08S*«nirS. 4) 
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CMJ. USWIKVMMTRIX HI .r.FTILI . IP. MI. mi, 4) 

pause 20 

CALL USUFHIVHHATRIX 01 *?.AL.Art22*10*miOS*NO. 4> 

LML IISURKIOHMATRIX TDi » 10*AHAT»10»M0M.«I.4> 
c ******* CONSTRUCT TOTAL STS! EH HATRIX Afll.DA 800000880 

MS-KCBStNO 

CAU. VttULFM8LA«*FLAn*W8*»fI**«*10*lO.AHAT*lO*lLRl 
WJ 718 1-1*03 
HO 718 J»1*IIS 

A! IAT U » J) -AHAT < 1 • J) ♦ ALAN C I • J) 

718 continue 

C CALI. EUiIcf<AHAr.l«*10,2.UEXO*Z#10*UKAR£A»IEJI> 

C CALL ll3UCV<]0KFl0F0VAL(iE(t0iUrX0fltS«X«4) 

CALL vm«.fK(W.AM*Fm2»WS*KX»ll0R8rJ0«10*FTlLI*ll>#IFRl 

DO 71? I-i »M3 

DO 71? J~1»MU83 

AliA f ( 1 • JlffS > "FT XL 1 C I * J) 

7iv umjtiwje 

DO 770 I-I»wm3 
DO 720 J-1.N0D3 
AHA T < I f « JtKS > -E083 C I • J> 

72ft .CUNT 1WJC 

DO 717 1-1 *110*3 
DO 717 J-1»N3 
AHATCIHIS* J>«FLUAf<©> 

717 corirzti:* 

IRUW*N3tN0&S 
CHART !>-• ATlt.DA TO* 
imo-'iM * 

C CALL (fSUFN(CHARt 13.AHAT* 10* IRON* IR0U#4I 
C ( AIL F10RFCAIMr*IR(Jtffl0»2«IICXa»2*10»llKARCA«IFR| 

C CAl.L IJSUCVC 10HEZ8CNVAL0E* 10* Wf ZQ»HS* 1*4) 

CIMWmi CONSTRUCT XFORH TO OET XIV COORDINATES til 
C 

C FIRST MODIFY L TO DECUHE T-L-LI13 
C 

IXi 721 I -1 .NOBS 
DO 721 J»1 *N0 


IF < J.UT . MSMiO TO 721 
IF « J.FO. X )L< I . JMffll •* L1IA1 < 1 ) 
771 COMTINUT 

c 

C Mill CONSTRUCT M INVERSE 
C 

DO 722 I-JUMJU 
DO 722 J~NStl*IkOW 
move i* j> -float co) 

IF < 1 . F(l« J) firNV<I*J>-FLOAT<l> 
727 CtXJTIMUE 

Wl 72* 1-1 1 RODS 
DO 724 J-I.NK 
IIINvaiNfi, 

724 continue 

c 

C CONSTRUCt n 


\ 


.-•Aktfk' 




DO 725 XMtIROM 

DO 725 .NNS41 » I ROW 

HI l *J> -FLOAT CO I 

IM I .tQ,J>mi..J>-FLOA7U> 

775 CONTINUE 

DO 724 I"1«RHDS 
DO 774 >*t *NS 

724. CONTINUE 

CALL WRJUTIl •N.NUKiNliRSf IO»10»FTai »10iXCX) 
DO 727 I-ltMODS 
DO 72? J"X*NS 
H(I4NS»J)-FT1L10*J> 

72? CONTINUE 


PAUSE 21 

C CAl.L (tSUFH(2HMt *2*M*f OtlROU* IfcOV«4l* 

C CALL 0SUFH(Sf8tIKVt»3*NtN0»S0>XR0U«IMW»4> 

cmmtt XFORH AHAT BY KFM H 888808808888880 

CALL VK&FF<HI»V*ANAr*lRGM«XROM»IRUU*lfr*10*FTILl*S0*IER> 
CAl.L V98S.FT (FT1L1 *H» 1R0M. IROMt IRON* 10* 10. AHAT* 10* XER) 


ettmm create. nem a -matrix for total, sylioi tmttttt 

«%AP<32.R£f>4)C<BCXItlJ>*IJ»l»Hl>*lI-l*N$> 

DO 720 1-1 •RODS 
DO 728 J-l.NS 

728 CONTINUE 

CALL VH!U ML • M.AH*NQ*8»lfS. MI * 10* 10*F TXL1 1 10* ICR) 

DO 72? X-I .MODS 

DO 72? -LM.NJ 

DC Xf RS*.I)-F1 IL1(I* J) 

72? CONTINUE 

CO 888808 CONSTRUCT NEM C MATRIX FOR TOTAL. 8YS1CM 000000000* 
RFAiH3?»R£C-3)<<C<II*IJ>*X>'fttN8)tIIM»M0) 

DO 730 I-l»ttQ 
DO 730 J-ltMS*IROM 
CC I* J) -FLOAT <01 
730 CONTINUE 
PAUSE 22 

CAl.L tiSWMOHAHATS *3*AMAT*10* IRUH»IR0II*4) 

CAl.L EIORFIAHAT* XR0M*X0*2fMCia*Z*10«MKARCA*ICR> 

CAII USULV< lOHFIUFNVAl.UFt lO.MEIS* IRUII* t *4) 

CAl.L 02^CmtOHFI»0«VCCTR»10*Z*IO«iROtf*IRtW*4) 

COOOOOOOO STORF HEM SYSTEM AHATtB.C.tt.Nf .NO 0(000000000000 
MR I TCI 20* RtC- 1 1 XROM. MI «M0 * X DOT » Zt RO 
URITI (20.REOJ) < (Bill* IJ)*IJ-1 .HJI* 11*1 clftOU) 

MR ( TL* ( :*OrRFf.-4 > i tCt 1 1 < I J) * I J- 1 • fRfftf I • f T- 1 • RO ) 

Itlf-IR(2U«20»1 

OFOIC (8IX T- 1 UT • F XI 'CMRRMT ' » ACCFSC^ ' DIRFCT ' »RTa - 1 02 ) 

MR I U.< XU I »RVO!i> UAIIAU II • I Jl • IJ*1 • I POM) « II-l * IMIMI 
COOOOOOOO FILTER INI UAL CONDITIONS AND XFORH TO FIND M(OI 0000000 
MRI1LC4* 731 ) 

73S FORMAT C1X. 43H CUTER INITIAL. CONDITION*. FIIR OCTOIN/A STATES) 
RtAil<3*e>(V0<Ii 1 >iI*1’NS> 

DO 832 f-l*Wf)BS 
Hi 832 

FTIL1U.JI— MI«vmMS*J) 





» 


*vr tN. 


\ 

\ 


S ' 




1 * •* • ' ' 7'W 


CM I. VHU.Ff (FTll I *V0* NOfcS*NS*t* 10*|0*V» 10* 1ER» 

DO 732 I-irNODS 
VO(IINS.i)-W<IIMCItl>-V<Itn 
AJ2 CONTINUE 

*33 ^{JwrtjlKtSlH USE THE FOtLOMING INITIAL CONDITIONS IN TINT REST 
CVUL </SUrN<5MX(0M*3*VO*lO*IIMJy*i*4) 

CtttIHtt* WJ Tint RESPONSE $$$$$$$$$$$$$*•$$$$$*$$$•*• 

incIlKAHLI) 

4AJ FWMArClK.2tiHWISM TO DO 1 INC REGFOMStlT I 

RKAPClrtlKK 

IMKK.I K.OHUI TO 444 
CAM. NOMC4 
AAA UKI1E(A»/.S4> 

7£A FORMA r<IXr40HUl8H TO REASSIGN OBSERVER AND TRY AGAIN?! 
READ(3r»>KK 
IFCKK.GT.OIUU TO 699 
C$$$$«$ RE6TORE INITIAL 8Y8TEN INTO 
RFAD(32«RFC*2)KS>MI»NO 

READ<;i2rREf;-3><<A<J*.f)..Hl*NS>*I-t.NS> 

READ<;S2tRf.l>4><<»<I.J»* 3-l*NII*IM*NSI 
RKA!»<32»RCf;-3><<C<I..II.J"I.»(S>tI-I»NO> 

URlTF<2<»tRFI>* >NS*NI *NU» IWITtZLRO 

yRirK<20rRE£.»2M(A<I*JH.Hl.«SnI-l»NS) 

U«nE(;!0,REr.-3HCB<I»J)»J-I.Nn#I-I»NSI 
UK t TE< 20 • Rtf>4 > < <C< I * J! * 3- I * NS 1 * 1»1 *N0> 

RE CORN 

ENQ . 


citmtmMiimitttMitmttiMitmutmiitiimimiiiMMMM 

SUBROUTINE TRANS 1 (A* IN* IN) 

RE At. A< 10* 10) t AT t 10* 10) 

DO 10 1*1* IN 
IK» !0 J'ltIK 
AT ( Jr 1 )*A( I* J) 

ID CONTINUE 
. 00 20 I«1*IN 

DO 20 J-lrlfl 
A<Ir.»>-Af<Ir J) 

20 CONTINUE 
' ^RFTUkN 
two 

CMttmmimiimmmmtmiMMtinimmiimmmmiiti 

cimommmNttmmiimititiitMmminTtimtntimuti 

SURNIHU INF VACT 
INfEGI R ORITJO 

RIAL VO<10»i0).LRrtRO<lO>»LIOWiClO>*VrSClO*t0l*V»CS<|Ofl0> 
lrAtm»lUOrIO)rRELaOrIO)*AIP«ARUO.iO>rEI(JOIf IlOtlOIt 
2VI2< 10* 10)*BLAN2< 10* 10)*R< 10* l0>*Wf?<lO*i)*IHM 10 1 10 1 
3.«tao*ioi»vomo f n.vF8T«puo. Hrvn:iPU<»»i)tAiAH22<io*io) 

ri'ai w'io.if»>*v<to* io>»viw<iOfi<»>»rao*i<*) . aiiamio.ioi 

RFAl. XX<1O*&0 )*Vh( 2O) »E(20)*X<2O) »I.RC< I0> iLlrit I0> »UJ( iO> 

REAL DI(j<20»20) *BI0IMV(20*20) iUKP((»< A*»0) 

RITA I. A(10»|(»)*B<|(»flO)»{;<IOr l<f) »Ut.»tRl!A( 1301 

COHNON/SrS/A*R*C»ZERO» I DOT* NS * HI *N0 
COMMON/ AIK3/F*AJIA I /F.IG/LlcE*LIM/|AR/Al /IJR/6 
COHMON/VEC/VAtE* X.UJ.U# XX* VtVINV 


CONHON/RO/OR flio* vo . VFS . VBCS » BLAN2 .LkORO • L I HKO . AL.AH22 

I2~0kIU0-NS 

IIHO 

I lii 'MItl 

I r < ARS (L IN ( xt) U .6T • ETROI ICHPLX* I 

cmmitumm create lambda-o i mtimmittiim 

DO 10 1-1*12 
1*0 io >i»i? 

RF.M 1 *J> -FLOAT <0> 

II<I.f(|..l) REM I *J)-LRC( IP) 

IO CONTINUE 

C CALL INiWEN* JHREL *3#REL* 10* 12* 12*4 > 

cmimtitmm create land da - 2 imumiiitmimw 

DO 20 M.I2 
DO 20 .HI *12 
ALI v HARCX*.l)>'AlAM22U*3» 

20 CONTINUE 

C CALL OSUFN(AHAI.PHAR*A»ALF«AR*IO»I?*I2*4> 
ctMtmmamn rel-aifhar tiMmmtmtmiu 
DO 40 1*1 *12 
DO 40 J*1*I2 

EIGO IF < I * J ) -KELT 1 * 3 )-Al FHARI I » J> 

40 CONTINUE 


C LM.L l*fiUFM(6HEIODIF*4*EIBPIF*l0» 12* 12* 4 > 

C PAUSE' 40' 

IF ( )( KPLX*NE»1 )G0 TO 34 

cmummii create limi mnmmittumi 


DO /I 1-1*72 
DO 71 3*1*12 
IHt.<I*J>”FLOATC0> 

IF< I .F0..1) IMMXr J)-CIH<III) 

71 CONTINUE 

CtlMimmil CREATE 116 IIMIIMtllllttilttill 

DO 72 1-1*12 

DO 72 .HI *12 

010(7 * JI-EI0DIFI I* Jl 

BICC1TI2.J) — Itfl(ItJ) 

»lQ(I*J4I2HIi2.(I«J> 

B10< 1412* JfI2)«CIBBIF<I*J> 

72 CON f TRUE 
122-2$ 12 

CALL LXHV2F<BIti*X22*20*BX6XNVtXB&l t a«K»XtttlCR> 
mi 73 1-1*12 
DO 73 3-1*12 
ALPHARt I*.f> n BXtiIHV< I» J) 

IHL<X»*’)«B7I)KHV<I*34I2> 

74 CONTINUE 
UO TO .IS 

Ctwtmill TAKE INVERSE OF EIGPIF $$$$$$$$$$*$$ 

34 CONTINUE 

CALL t.INV2F(EHil*IF • 12* 10.MPHAR* IDGTrtR.AREA* ICfcl 
£$$$«$#**$$ FREMILT BY V12 •$$$$$$$$$$$$$$$$$$$$$$$ 

33 CONTINUT 

DO 41 XM*N8 
no 41 3»1«I2 
V12(t».n-VO(I*JTNS> 

41 CONTINUE 
C r*AUSF/41' 


ro 

o 


ORIGINAL PAG” IS 
OF POOR QUALITY 





c iMix usuFfMTHiNVERSE* /fTtiPUAR* 10*12* 12*4 > 

CAI L VHULFHV^'AU>|fAftrM»>I2tI2*10»i0*f*t:t»lC’rI£fCl 
C PAUSE 'RKL' 

CUMUUM POSTMUT BY BLAM2 tltttttmttmtllt 

C-ALL WttJLFMRELtW. AH?*HS. I2*N1 . 10* lO* AI.I1IAR* 10* TEN) 
t PAUSE 'M AM2' 

FOS1HIN r BY F ttMIttlMMIttlttMl 

CALL VMMI.FF<M FHoR*F*NG*NT ,NS» 10* 10*REL> 10* IER) 

C PAUSE 'F' 

RESULT TO VIA ttllttUttlMlItt 

Ml 50 I«1*NS 
DO 50 J*1*MS 
R< i * Jt pjt 4V»< r > j> 


50 CONTINUE 

CttIKMMt FUST HULT1PI.V BY ASSUME.* EI8CNVECT0R lltltttlt* 

DU 60 1*1 *Hf# 

VGft<I*i)«U<I»I(l) 

60 CONTINUE 

C CAM (jSVFM(iHVQk*3*VUR*l(»fNS»l«4> 

C CM L IISUFN<SHR*i*R» tO«N9»NS*4l 

CALL VNULFKR*WR*N6*NR*i*10*10*VFIi1HP*lO*IER) 

IHI 70 1*1 *KS 
VFS <1 * IU I -VFSTftFt I * n 
70 CON 1 1 HUE 

lFUCNPLX.faM) SCI TO 111 

ctttttttt*** prfmult by vi 2 

CAl.L VMULFKVJ2* I HI. » NS* 12*12*10' 10*REt *10* IER) 

cmmmttu* pobthui.t by blam 2 s*e«**e**«tt««t*«* 

FAIL v«m FF<R£L*BLAM2*NS«I2*N1 *10*10. ALFNAI* 10t XER) 

postkult by f 

CM. I. VHU1 FF< At PHAI . F» NS. Nl»NS* 10* 10. Ml *10* IER) 

cmmmtmt compute eioenvectors 99999999999999 

DU HO 1*1 *NS 

voi<r*i)*v<i»ia4S> 

HO CONTINUE 

FAIL VttDLFC<HI*U<II»NStNS»l*10*10*VF8TMP*lOr IER) 

DU VO !*l*Nft 

Wt>(I * X0)*VF6<I*TU)-VFBTHP<I»I> 

»0 CONTINUE 

CALL VMUl FF<R.VOI.NB*NS*l*i0*10.VFBT«F*10*ICR> 

FAI L VMU1.FT <HI* Wft»N£»NS* f *10* 10* VF81P* 10* XER) 

DU 110 1*1 »M8 

VFSC I » lUf 1 )*VF8TNP< I • K >4VF8TP< I * 1 ) 

110 CONTINUE 
1Q-104J 

111 CONTINUE 
lCMPl.X-0 

IF < JU.LT.NS) SO TO 1 

RF TURN 

END 

CMMtmmMtMimimmiMimtmMftmttfmiimimmi 

ntttttmtmmmtMmtmtMiitimmmimmiimmim) 

KUPRmiT I HE RULTJST <C J) 

INTEGER (IRISH 

REM. XX< 10*10) •VA<: , 0)*E<201*X(70)*LKt:<10)*LZH«f0)*tfJ(t0) 
REAL U< 10* 10) *V< 10* 10) *VINV< 10* 10) *AI.< 10* 10) *AI.AH22< 10*101 
REM. VKft< 10*10) *VDES( 10* lO) . BLAN2I f O* 10) » VO < 10*10) 


RfAL AUO*lO)*B(lO«iO)*C(l(i*l<rT*LR(»Q(i»;n0rctillF) 
i:0HNnH/SY8/A*B*C*ZER0t ID6T»NS*NI •MO 
i:OHNOM/VEC/VA»E*X*tfJ»H*XX*V>VIflV 

COmOfl/RU/Tmimi* VC) f VFS e VDES *BL AH2* LRtlRU»UURG *AL AH22 
LONMON/E Hi/LREtL IH/PAR/AL 


citmmt calculate AcruAi v in fuil system mmimiMmiiu 
CAM. VACT 

C FAI L USUrN<3llVFR*3*VFb*10*NS*NS*4) 

C FAIL USlfFNMHVDESi 4*VD£S» A0*NSiNS*4) 

CMttimi CALCULATE COST FUNCTION tuimMimtiMmiiimint 

CJ*0 

J*0 

3 J-J41 

1-0 

10 1*241 

C J T F HF« < VF8 FI * J > - VDCG < I * J ) ) 1 1 2 
C J ! F.) If MP4 AL < I » J ) 4C J 
1F<ABS<LIM< J) ) .LE./EKO) 60 TO 20 
CJ-LJ4LJTFHFBAL <!•.:» 

C J-C J4 < VFS< 1 1 J4 1 ) - VBES< 1 * J4 II I t«2*2*AL < I * J4 1 t 
IFCX.NE.NS) 60 TO 10 
.1" J 9 f 

20 IFU.NE.HS) 60 TO 10 
. IF<J.NE,NS) 60 TO 3 
RETURN 
. END 


cmtimumiimmmiiiiiiiimttmmimttiititimiiititit: 

SUWHKIMNE R OUR AD 
HJTtOER mUUO 

REAL AUX'1< 10*10) *MIX2< 10* 10) »MIL4< 10*10) *FV(?RE< 10*1) »FV0I< 10*1 
1VAUX1 <10*1 ) * VAUX2( 10* 1> *B£fAR< 10*1 ) *?E TaRC IP* 1 > * 

2VU <J0*10>*BErAI<10*n.2ETAX<10*l)*RR<10.1) *00(10*1 > 

Jekf) (10*1) *0ft( 10* 1 )«()< 10*10) *AL< 10* 10> 

RLAL VG< 10* 10) * LRURU( 10) * LIUR()( 10) *VFS(10«10> *VDE!» ( I0* I0> 

I * AL41IA1 < 10* 10) (RELUP* 10) * ALFNAR< 10* 10) *EZUDir( 10* 10) • 

2V12< 10* 10) *BLAM2(10f lO) »R( 10* 10) • VUR( 10* 1 ) * 1ML < 10* 10) 

3»NX< 10* 10) « VOX < 10* f ) « VFSTNP< 10* 1 ) r VFSTP< 10* 1 ) 

RCAI MI (10*10) *m.<10i ID) *N1C< 10*20) *FLC< 10*20) iNLC<10*2O> 

RE At STAR <20* 20) *M.< 10* 20) *81.(20*20) 

REAL V(1O*1O>*V<10< 10)* VIMV(10* 10 )*F( 10*10) • AJIAT< 10* 10> 

REAL XXT10* 10) * VA(2O)*EC20) *X<20> *LRE< 10) *LTH< 10) *U.I< 10) 

REAL A<i0*10>*B<10.10)*CU0rl0)*Ut»/tREA<lJ0) 

REAL 010(20*20) * BIO INV< 20 *20) *MLBIU<460) »Al.AM«2(t0t 10) 

LOMMUN/S YS /A » B . K * ZERO * 1 1<0 T • ITS . II X * NO 
* C0Mf1U»/NSFA/ftt.*NL*«I.C;PLCr«lLrSTAfe*flL*RL * * * 

COMMON/ MIX/ AUX i * MIX2 »AUX3 
LONNI itt/MM/F * ANA T/E lU/LftT «L1 N/PAR/AL /uR/O 
COMMON/ VELYVA * E , X * IM * U * XX » V. VI KV 

FOHMOM/RO/OR 1 01) * VO * VFS . VDES . W. Att2 * LRORG , L X ORO » ALAM22 

I2-URIFU-NS 

1FLA6-0 

IK) 10O 1*1 *NS 

DO 10O J-'lrKl 

C(lt J)^FLOAT(O) 

100 CONTINUE 
I I-O 

110 Il-llfl 




ufiismc 



» * 



IQ-0 

i?o iu-um 

III* 1(2120 

iriABsaimmn.cir.zrRoiao to 34 

Rl'ADC III. REC-3) ( (ML ( I> J>* J-l >NI ) *1— 1 »K3) 
KCAnmi*REc-4>(oiLfXtJ>*j-i*fji>*i-i*Ki> 

IH» TCI 35 
34 IS-NSTNI 
NI2-7*W1 
NS2-NS*2 
IMS-WSM 

hl*ad< ui.Rtx-3) t meat j>* j-i * is>* 1-1 *ns> 
Rf.A|UIU>REC«4M(HCUtJ)»>Ii29I(IMiNS) 
RrAnauiRt‘c«5M(ntc(i»j>tj*i»z8)fi-i»Kn 
W AH< III'REOAI < cm. ( 1 . J> » J-l *NI2> • 1-1 *NS > 

R1 AIK lUf REC-/M <w. C I* J> $ J-l *MI2> 1 1-1 # MS) 

jts continue 

C PAUSE 'PFX' 

CALI. PEX( IlrlUrlELAO) 

C PAUSE 'PFX' 

CALL PVPdltltt) 

C lAUGt'PVr' 

ICC ADSCLltK 1(1) I •(.£ « ZEROIGO TO 125 
IH.A02-1 
IFLAHI-I 
125 COM T INC C 

CtltMttmmt CAtCULATE ALPHAR ********************************** 

uiimmmutt create iahrda-o i titimmmmm 

W» 10 1-1*12 
DO 10 J-i *12 
RCl.< I* JI-I'I.IIATCOI 


IMI.EP.J) REMI»J>-LRE(IU> 

to continue 

CREATE LAMBDA -2 ******* **************** 

DU 20 1-1*1? 

DO 20 J-l *12 

ALHIAR C X * J ) -AI.AH2? C I « J> 

20 CON I INC IE 

C**************** REL-ALPHAR ********************* 

no 40 i*.i*i2 
no 40 j* 1 » 12 

ClUDIK<(*.l>-RLL<I*J>~AtPHAftCIf J> 

40 CON f INC Hi 

ir liruuu.MM > lo ro 134 

c************ one a ft. i.ih*x ***************** 

DO >1 1-1.12 

no /j j-i*x 2 

»MI. C I * JI-PI.CIATCO) 
irci.Pii. jiiHt.cif Ji-i.iiidii) 

/I CONTINUE 

C**t*t*******«tf CREATE BIO ***************** 

DO 7? 1-1*12 
DO 72 1-1*12 
|U0C1*J»»ET01UFCX«J) 

Dio(ifX2*j)— ira.<i».A 
IUGC1»JIX2)-XIU <I*JI 
. D1CCVII2. JU2>-EX0DirUf J) 



72 CONTINUE 
122 - 2*12 

CALL LZNV2M»ia*I22»2OtilQIKVtIMTtllK»I0*XBt> 

DO 73 1-1*12 
DO 73 3*1.12 
ALPKAfcl I » JI-BXQINVC I * J> 

73 (CNN NUE 
CO TO 134 

C*********** TAKE INVERSE OF EIQDIF ************* 

134 CONTINUE 

CALL LIMV2E<EI0Dir *12*10. ALPMAR*1WJ1 .IWCAREA* IER1 
C********** PR1MLT BY V12 ************************ 

134 CONTINUE 

DO 41 1-1 *MS 
DO 41 .1-1*12 
V121I* JI-V01I. JtNS> 

41 CONTINUE 

PAIL VNULFF<Vl2*Ali > «ARfNS*12*I2»10*10»lca*10*IERI 

c********** posrnui.T or w.wu ******************* 

CAI L VKI» FFCREL.W.AH2*KS*I2*W1*10*10.AI PHAR*10*IER> 

c*tm««mm calculate re tar ********************************* 

DO 130 I— 1 .NO 
rVUREH* 1 >-AtOCl C 1*1(0 

I30 continue 

CALI. VI(ULFF<F* PVPftE *IH • «S* 1 * 10* 10* VAIJXI .10* IERI 
DO 135 1-1 *KS 
VOKU *i)«V<X*I()> 

135 CONTINUE 

CALL WIN. FF<U.WR*Klf NS* 1* 10* 10* VAIIX2* 10* IER) 

DO 140 I-l.NS 

BLTARl X . 1 >**VAtIXi <1 • 1 MVAUX2C 1 . 1 > 

140 CONTINUE 

C************* CALCULATE ZETAR ********************************* 

DO 143 I-I *WS 
DO 145 J-lcNS 

vn<i*j)«voci.j> 

145 PON fill* IE 


«*** ************* cmS 

.so 

imFLMII.NE.O Ml U) 1*7 

£!2IJf!**!****** -i.'tiAfE au'iiai imtmiinuitwMiiititin 

rnciHN.r by vis ****»t»t*tttt***tt*t*ttt 

cul vwii.iT<via,nu .ns.ia.iz,io.io.REt.ic..icii> 

citmmimii iosikult by blam 2 *tt*tt***t«**tt*** 

I'M.i. VHM.IT < Rtl. , BI./lMi! . HS . 1 1! . N 1 . 1 0 . 1 0 . ALH l/t l , 1 0 . 1 1- K ) 

etilittititsfittit ifjnnuiK wtiai 
dii 

PVCIX (1*1 1-AUXK ItlUf J > 

155 CONf /Nl NT 

PALI. VMlW.FF<F *PVPX tKl iNS*l»lO. JO.VAIJXi . 10* IER) 

Dll 1A0 1-1 *NS 
VUICI»l)-0(f*I(j|f > 



rr^iniv; 1 




I AO CUNT I HUE 

CALL VrtUCrECU,WI,NI#NS*l#10rl0.VAUX2*l0,l£R> 

Ml i AS X-I.NS 

&t T A I ( X , 1 l-VAUXJ C I * 1 > FVAUX2U . I) 

IAS CONTINUE 

ctttttmmti compute zetai ttttttMt****** mmttmu 
i:am vmijuti vi t .pvoi tM3#NS» i # 10 . 10 # zetai « io» ieri 
cmttmtmtti calculate uu ittummittMittxtumiM 

CALL VHUtFFI ALPIIAI tBETAX »NSiNI»1f 10» 10>RO» lOr XCK) 

calculate tw 

CALI. VMUUT <ALPHAR» BETA I »NSf HI r X. 10. 10»RR» 10# XER) 

WJ I/O I-l.NS 

URUilWETAHItlHORUtl) 

I/O CONTINUE 

c ttimitlt BEE XF THIS IS Sf.CUMP TIME THROUGH FUR COMPLEX t«l 
ii ori.Aoz.KC.n ou to iao 
I A/ JJ-III 

GO TO 1AV 
1A0 JJ-IUI1 
1AV CONTINUE 

cimmmi cuhputc gradient for real case 
ip-o 

l/S IP-IPM 

G< I ! t.JJ)*GC XI r JJ) M0FS( IF* IRI-VDEUI IP* 101 It <RR< IF# t I 
lt2tAL < IP* III) I 
IFOP.NK.NS) GO Til 1 73 
IFIIFI Alll.NE.tIGtt TO 203 
CKIIttmt COMPUTE ORAM ENT FOR COMPLEX CASE 
IP-0 

too ip- in i 

Gill » J.f )-0( II» JJ) - <UF8(XP»Itf)~ VDEHI IP»XR) ItlfrtH IP* 1 1 
) *2*Al. < IPr lO) Ml <VFG< IP* X(Hi ) -VDESIIP* lflti ) ItlflR? IP* II 
211101 IP* I ) )*2tAL< IPt IDH > > 
irUP.NK.NSI GO TO 100 

cmmtm check if this is first time through immitmmi 
ir< in AK2.KE.il GO 10 201 
If l AG2-0 
Ml IU5 I-I.MS 
30 I OS .1*1 *NS 
tUIX J < I * .|>«AUX2< X t .11 
105 CONTINUE 

1*0 1V0 X«i*NX 
DO I VO .fait MS 



O O 

-n g 

-o o 

o '■? 
o > 
w f- 


iC *cj 

c: ');> 

> 

r 


1 


1V0 continue 

GO TO 125 * : 

201 CONTINUE ^ 

IU-IUI1 ! -S 

203 CONTINUE 

IF< I0.C f .NS) GO TO 120 
XI (II .I T .NX) GO fll 110 

CALI. GFNORM<NI«NSI 4 

RETURN l* 

LNU ' |j 
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, ************t****t**tt*t*t*t****tt**t**t****t***t***********t*t******* 
********************* SP+rCTRAI. ASSXHNKKNT PACKACK ******************** 

ENrt:r< »esxre» mine ok cperat i ok > H oritr^x* . i » a » , . . , vi 

r t 

************************** none iiuata entry ************************* 

**********ENTER OR CHANCE SYSTEM PAR AKE IE RS S 


PREVIOUS VALUES? 

? X 

N '*“ ' ^ NO" 4 1 NOT- H ZERO** .OOOOOO+OOOOO 


WISH TO CHANCE? 
? 0 

MATRIX A » 



1 

6 

a 

7 

A 

4 

s 

t 

I54000E+00 

--. 320000 c oi 

. 1S4000EI01 
• .OOOOOOE+OO 

-•. 420000 E -02 

•OOOOOOE+OO 

-.744000+: +00 

a 

-.VVAOOOE+OO 

•OOOOOOCIOO 

-. 117000ET00 
.OOOOOOE+OO 

-.avsoooE -o;j 

.;mooo+:-oi 

.200000E-0I 

3 

1 

.24vooor+oo 
-.1 J2000CI0] 

-.saooooe+oi 

.OOOOOOE+OO 

- . XOOOOOE+OI 

.OOOOOOE+OO 

.337000E+00 

| 4 

j' . 

.000000+: 100 
•OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

. iOOOOOE+Oi 

.OOOOOOE+OO 

.OOOOOOE+OO 


.OOOOOuE+OO 

.OOOOOOCIOO 

.'ooooooe +oo 
.OOOOOOE+OO 

.OOOOOOE+OO 

•OOOOOOE+OO 

-.aoooooE+oa 

6 

.OOOOOOCIOO 
- . asoooot +oa 

.OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

.OOOOOOE+OO 

7 

WISH 

. SOOOOOE 1 00 
.OOOOOOCIOO 
re CHANCE? 

•OOOOOOE+OO 
-.SOOOOOE +00 

.OOOOOOE+OO 

.OOOOOOE+OO 

•OOOOOOE+OO 





ORl&NAi: n 

OF pool? Q ; JAUl« 

* 

MATRIX 

» t 

I 

*» 



I 

•OOOOOOEfOO 

•OOOOOOEIOO 



a 

•OOOOOOEfOO 

•OOOOOOEfOO 


* 

a 

•OOOOOOEIOO 

•OOOOOOEfOO 



4 

•OOOOOOEfOO 

•OOOOOOEfOO 



S 

•aoooooEfoa 

•OOOOOOEfOO 



A 

• OOOOOOL f 00 

• asooooF.foa 



7 • <><>00001' 1 00 

WISH TO CHANGE? 

? 0 

MATRIX C J 

I 

A 

•OOOOOOEfOO 

a 

7 

3 A 

4 

• 

!3 

i 

. loooooEfoi 

•OOOOOOEIOO 

•OOOOOOEfOO 

-.lOOOOOEfOl 

• OOOOOOEfOO . OOOOOOE f 00 

• • OOOOOOKf <>(> 

a 

•OOOOOOEfOO 

•OOOOOOEfOO 

•OOOOOOEfOO 

•OOOOOOEfOO 

• lOOOOOEfOl •OOOOOOEfOO 

• OOOOOOtifOO 

a 

. OOOOOOEf 00 
•OOOOOOEfOO 

. lOOOOOEfOl 
•OOOOOOEfOO 

. OOOOOOEfOO , OOOOOOEfOO 

• OOOOOOE 4- 00 
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4 .OOOOOOEfOO .OOOOOOEfOO 

« 0000 JOE H>0 .OOOOOOEfOO 

WISH TCI CHANGE? 

? 0 

wish to cxrr krum this mode? 
f i 

*****************At*A*****t EXITING MODE 1 
TERM mate this run or select next mode: 

WISH TO TERMINATE? 


000 00 OUT 00 


•xooooocfoi 


•OOOOOOEfOO 


************************* 


/ ****************************************************** 1 ):*************** 
********************* SPECTRAL assignment package ******************** 

ENTER DESIRED MODE OP OPERATION. M0DE=0. i , 2. . . . , 91 
T V 


WANT TO ENTER NEU 
0 

V MATRIX! 

1 

A 

GRIG I NAl EIGENVECTORS ? 

2 3 

7 

4 

S 

X 

•aaaoaiKfoo 
. 1 2V240E-02 

uvooaii+oo 

•OOOOOOEfOO 

-.Bsivoai: -<>i 

• a74401Ef 00 

•37417AE-01 

;» 

.XUA40Ef00 

• ssitioai:-o 4 

•30VA71Ef00 

•OOOOOOEfOO 

• Aasvoyi: -oi 

. asAasoE -oi 

.ftA444!vK-'03 

;* 

-•.loaaaaEfOi 
• 46o/.sat: -ot 

-.asvvsacfoo 

•OOOOOOEIOO 

•AOSAyacioi 

- .Aas.viai: 01 

-.i7yyoii‘-oi 

A 

-.aasyovi-foo 

J06AA1E- 02 

•02VA27KF00 

•OOOOOOEfOO 

-.3737S iCIOl 

• 7ia;5S7Ef01 

.«yv«oai:-o:j 




5 

,0000004.400 

,0000004.400 

.0000001- (00 

.0000004.400 

. 1000004.(01 



•OOOOOOEIOO 

.0000004:400 





i, 

• O0000OEI 00 

.0000004: 400 

.0O0000CIOO 

.0000004:400 

. OOOOOOEIOO 



,i0O00oi: hu 

.0000001:400 

» 



J 

7 

. X8AVV7E “02 

X 502t51E 400 

.Z2755AE-01 

• 27V404C (00 

-.V5V42AE-03 


• 

.7437S5SE-04 

.1000004:401 




j 

wish n 

CIIARt.E ANY 

VALUES 04* VY 




? 

? <• 






i 

alami 






j 


1 

2 

;» 

4 

a 



6 

7 





1 

• 881901 E--01 

.12AV48E401 

.00OO00C 400 

. OOOOOOEIOO 

•OOOOOOEIOO 

i 


,00000011 (00 

.000O00C (00 




J 

j 

2 

•12694UE401 

-.881901K-01 

.0000004:400 

.0000004:400 

•OOOOOOEIOO 

1 


• 000000C (00 

.0000004:400 

* 



1 

;» 

.0000001:400 

.0000004:400 

-. 10854SC4O1 

.00000045400 

.OOOOOOEIOO 

i 


.00000OC (00 

.0000004:400 




1 

4 

,00000011400 

.0000004*400 

. 000000C 4 0O 

-< V14482E-02 

•OOOOOOEIOO 

t 


,O00O00C 400 

.0000004:400 





s 

•OOOOOOEiOO 

.0000004:400 

,0000004:400 

. O000O0C400 

- . 2OOOO0EI02 

; 


,OO0O00C400 

.0O0000C 400 




• ♦ 


.OOOOOOEiOO 

.00000045400 

.00000045400 

•OOOOOOEIOO 

•OOOOOOEIOO 


* 

•2SOOO0C4O2 

.OOO00OC4OO 




i 

7 

.0000004:400 

.0000004:400 

.000000C 4 00 

•OOOOOOEIOO 

•OOOOOOEIOO 



.0000001*400 

-.5000004:400 





WIOH 10 III SPLAY K(i 

HOPEI. T 




i 

T X 






-I'.: 

(vt.ljUCL.lj 

ORPER HOPEI 





♦ 

i 

4 

MAllUX 

a: 





* t 

t 

. 1 


1 

2 

J 

4 


1 

■ 1 

1 

•aoxyoiE-oi 

.12A94&C401 

•OOOOOOEIOO 

.OOOOOOEIOO 


K 

2 

.1209404:101 

-.881VO1E--01 

.OOOOOOEIOO 

.OOOOOOEIOO 



.» 

,0000004:400 

.0000004:400 

- • 1 08545C 4 O I 

•OOOOOOEIOO 


4 

4 

, OOOOOOEIOO 

.0000004:400 

•OOOOOOEIOO 

-.VX4432E-02 



MATRIX 

Oi 





j • 


X 

2 




\ ‘ 

I 

.XMU'JStiOX 

-.1444404: -01 




1 . 


,Aiova;5C40o 

. 722V0VE-01 





3 

•27/21 5E40O 

-.28780/4:400 




•j 

4 

• 24?75.(E(0O 

-.1527514:400 





MATRIX 

CJ 







1 

2 

:i 

4 


$ 

1 

.381 151 Cl O0 

.3124V1L* 01 

- . X57V44E TOO 

-.5125881*- 02 


‘ ; . 1 


I 




ORIGINAL Pi;» : ;j? 
OF POOR QU^'J n 




in 


~.io.j.12/kioi 
. U1A4UU00 

.22SVovkioo 


--..ISVVSiJKlOO 

..tova/ieioo 


.iosav.jeioj 
. a::svove oi 

-..I7.J73 1KIOI 


-.A»SAi;*C OI 
.2:5a.isiie-oi 
• iM3.<S.TI0l_ 


- — — — — ... »-i- ZkfOO *• ■ if tin »7iK<K *rnn 

Hlflti!;!!******* «” C*«OW*UIE A^nwi«inmM*S-*i 

mtltlUMIMMMtl MtiOE 2iF.IOCKVM.llK ASSIGNMENT t«Mt*3*4***M 


ENTER (111 CHAN Oil EIGENVALUES I 


MCVIOUi; VALUES? 

F 0 

LAMRDA it 

RLAL» < OOOOOOE ( 00 


1MA(I» 


Ml Sill r o chanoi: T 
f 1 

CNfi.it new VAu»:<8> t 

r ••*.« i.s 

LAMMiA .‘.'iKCAL* 1SOOOOKIOI 
NEXT EIGENVALUE! 
fkCVIUU'J VALUES? 

> o 

LAMUDA .11 

ItCAt* .0000000 00 I MAO' 


.OOOOOOKIOO 


•IKAS* - .ISOOOOCIOt 


.OOOOOOKIOO 


-.looooodoi 


Ml OH 10 CHANCE? 

» ( 

CN fl'R NKU VALUE <«} I 

r -2 i 

LAMUDA 4iKKAt.» - .i.’OOOOOKIOL • IMAG* 

MltH TO EXIT KROH THIS HOOF.? 

F 1 

F.NTI.r< DESIRE I) EAR f I Al. KIliENVKCflM A^I RCNf ******************* 
KIGLNVECTUR V I 
F 20 0 A J <1 <1 O 0 

eigenvector v j 

f o o 0 o :<i o -a -4 

UCK flic EMI.UUINC v MATRIX FOR INI f I Al. ASSIGNMENT 
RLMLMDEK Mill Ell V ARE COMPLEX 
INI HAl. GUESS EUR Vi 

» 2 


1 .4Hf.V0.JKI 02 

2 ■.IWS.'OllOi 

.1 . i 24I2IKI02 


«f.V!Js;ioi:ioi 
. 1V403/EI02 
•31V0ISKI01 


• .fl.JVAU.100 
-.VAVSiJ.IE Oi 
.70V22U: Oi 

mmumurnuu mm * r 

********** ENTER (JR CHANCE EIGENVECTORS! 


3 

*Ei0O 

“.vtioosunoo 

.47(l8Vr,EIO( 


.IliHOV/.EIOJ -.20//32KIV0 . llfli'tl/EKx 



TTtLVIUUO VALUES? 

F 0 

CICENVKCUilt V I! (REAL) 

.OOOOOOEIOO 
•OOOOOOKIOO 
.OOOOOOKIOO 
•OOOOOOKIOO 


(MAC) 
.ootioooi: 1 00 
. ooooooi; 1 00 

•OOOOOOKIOO 

•OOOOOOKIOO 


UISII f£1 CII/iNCKY 


«vi< mt> 



ORIGINAL PAGE IS 
OF POOR QUALITY 




t i 

cuter a mew desired vector t 

f 411./ 6.4*6 -1.3/ IV. 4 12.4 3.4V II. Itt -.20H 
cimrccx vdi 

.4H/O0OCIO2 .i.I/OOOEIOl .1240001*102 .HlHOOOE 101 

.1V4000EI02 .34V000EI01 -.20H000EI00 

ACTUAL VECTOR > 

.4H22/SEI02 . 12IOV 4EI01 .12/0161*102 .ftOOl.'IOEtOi 

• ivoi/.ieios* « 3/366 heioi «.<» 37:o/i: -02 

ERROR VECTOR 1 

.4?2b46E*00 -.2SH0V1EI01 -.J01S/3LI00 . t/fl/OlEICO 

-.4132X/EIC0 -.2466H1EI00 -• 1VV62UEI00 

i.kixi fit or the j;i::ureb vector » ss.osavki 

LniOTlI or Till* PROJECTED VECTOR- S1.V(L.'U4/ 

urnofM or rue ckrdr vector ■ 2 .H 34 VS 5 

10 fill’ ERROR ACCEPTARtEY 
t J 

ci hervci: tor v 21 irerd < in aid 

. 40 /OOOE I02 - . 6V6000E 1 0 1 

-. 137000EI01 -. 1V4000E to:; 

.124000EI02 -.34VOO0EIO1 

•H1R000EI01 . 20B000EI00 

or XT CXUENVECroiU 

tlOENVCLlOR V 31 (REAL) IIMAO) 

. OOOOOOE 100 . OOOOOOEIOO 

.OOOOOOEIOO .OOOOOOE 100 

. OOOOOOEIOO . OOOOOOE I 00 

.OOOOOOEIOO .OOOOOOL1 00 

VICK TU CHARM"! 

r i 

CMfl.R A HL'U DESIRED VECTOR I 

r -..106 .344 -.VIII -.09/ 4./V ,070V 1.4V -.3 
COMPLEX VO I 


.6V6000EI01 
, 7H6JVVEI01 
- .V03VV.lt* 1 00 


-• , 306000E 100 -.VOIOOOEIOO .4/9000EI0J . 14VOOOEIOI .344O00EI00 

-•.V/OOOOE-Ol . /OVOOOE 01 -.OOOOOOEIOO 

ACTUM VECTOR I 

-.414202CI00 -.7006041100 .4H.UA0EI0 1 .13VVS3E101 .3SOV1SLIOO 

• 203/3/C 1 00 , 12/41 VE 100 -.4V30//EI00 

ERROR VECfOR 1 

. 1002021:100 -.2»o.tv6Eioo -.434oivir-oi .voaaase-oi -.6ViS32K-o;* 

••• 3O0/37C 1 00 -.S4S1V4E-01 -. 6VJ320E-02 

I. ECO III or DIE DESIRED VECTOR S. 13/00/ 

I.Elllilll 01 TOE PROJECTED VECTOR" 3.13UV3X 

LERIS I II 01* TOE ERROR VECfOR ’ .44062:.* 

IS THE ERROR ACCEPTABLE? 

' 1 

ciuERvrcnm v 4t irealj iihaid 

- . 304O00EI0O - . 3 *4 OOOE IO0 

- . Vll A OOOE 100 .V/OOOOKOl 

•4/VOOOEIOl *./OVOOOE 01 

. 1 4VOOOEIO 1 .300000C 100 

««.«««. ...CORTERTS 01 * C! JRKI.'T • DAIA EILE INCLUDE: 

MATRIX V I 

1 2 .1 

4 


I 


. 4II22/4S4 1 0//2VE 1 02 
. .ISOViS.llSI 4102E 100 


. S'll43VV:i.UiS!1743i: 1 0 1 


UlC’OiVSJSVAOUTOO 



. l210V3ll.l2.IU3i 21" 10 1 


.xv«i 3 ::s./ 40 i;:s.hto 3 


0.40 S *60/ f 1 ivrisso 




. /o.t.MAAAWtAA/.'Kxoo 

.ja/HlS/.llK/X'OOKXO:! 

. iiV 4 j v.ivHssa-JstTXoo 


original PA'r: :r 

OF POOR QUALITY 

• .(.maaU'MS'/a/ 44 I.»ox 
- . HJ/aoAVA/VA IS/i: -o:: 


4 .H0<H/V8/0:’A/V::E4<!I 

• ■•4?3C74/V/Bt*10.IC 400 

u<:;m rti oiyi’i.AV rur mnn/tUfieo nitn.nvr.ur.Koy 
o 

OAIM MATRIX rt 

X 1 


< . 1.J4 V00tto0HX;:42t‘Hli 

i - . .WA.tO A.5HS.I 4AAVK 4 0 1 

,3t i 8VJ4HA4HOA/K40a 
MATRIX AHA f i 

1 


1 *.ai7/75iJ»H/;V/F40J 

,2 SO iOlMH:J2tJ»KiOt 

2 ~.742VJ02Oi-4OKV/KtOO 

.iaon 5 » 4 ; ,, /.»i 4 y;»KKu 

.1 . /4A X All/X SHAUSAK * 00 

- • XMVHA.ISO IMS'! 1,'KKil 

4 .2J*2Ui2OA.MX,’KtO0 

- . 4Jon.i:u vss/hvskko i 
wish ro exxr krom this mows' 


■.lASM.t/UJA 1.114. IOi 
• A0A.J4/O IVOH'.KIVK 100 

•* 

* 

. avsSMSo 1 .’va xacku 
X0S;U!S4420lliUVi:4OX 
. 2».trV043S50A/a0K *00 
. .VtlHS/HOSSAASOt: 400 


. 4»;WA0i B/a06.t4l'401 
. X .»VV8J3f.4SV f2 * r 40 X 

.1 

. x Ava.’/vvA ivs//x;4oo 
- . 8tlOH!5.IH//AO/ I Bit 40 X 
3 

- . X HOBOS 10/4 7/044! »00 

- . .ixAiHH'/ajs.u /si: too 

. S.W1S/7A 1V5X 1 AX' tOO 
.HI IBV/3//HBVVA1.400 


*«**«*»***»*»■** KX/Tlf.'H Knot* .1 *«t***********«********** 

MIN tUitktX X'l 

i :: ,* 4 

X .UlVOXKXOJ -• i AS14.1X* 401 . X AVa/HEXOO - . 1/2JU/K 4<U 

i AU6J06E4OI .A04.M7K400 -.SOOHWlKloX . ,t< XHV.IK 40..’ 


iw»tr rx ox's 



1 

2 

.1 

4 

4 

. XV4.<S4K4«a 

. a.i.1a.1::x'4o« 

XOaVOAKXOO 

-..iAOjaof-<n 

a 

•A/A4A0K40< 

./24/24X:40X 

./s/isvxr-oi 

.V/AOS/K-OX 

.t 

.XOAlS.'K+OI 

-.V/40-10K-01 

./Ol VH0K40:X 

. U.U.5/M00 

4 

WANT 

- . 4:.'4a.’SK4O0 . A.is l SSK 400 

ro <.'ONT<M/K SKARCHt 

-.HS4IH0f4Oi 

-.:ivx.i/Sfc*4<n 

«*»*t**«**t*t**t*t** 

WKIUHTIRH CONSTANTS 
wxrniHfsi 

RH KKiKNVKCrOR 

(M4'K<|VKHKNT******««**««t***t*«* 


1 

2 

.1 

4 

« 

. I0<I<IO<IK4<I < 

• <<MMMMIK4<H 

. <00<MIOK4<H 

. <<MMI00K4O< 


. <<MMI00K4<U 

. <00<MMIK4<H 

. <<MMMMIK<<H 

• < 000004:40 i 

.1 

• <O<MMI0K4<H 

. <0<I0<MIK40I 

•lOOOOOKXOt 

• I<MIO<MIKKI< 

4 

WISH 

. I0<MI<MIK4<H 
fO KHmNOK' , 

. <00<I00K4<H 

. i <MMMM<K40< 

. <<MMI00K4<H 



•*t T "HU llffl/ 



/ 

/ 




/ . 



•/ ’ 


9 
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ENTER NEU VALUES l 

f * i ins iK» i i its iks it/ »:/ n u:/ »r/ t a 

COST* ..MJVV.O+OH 

GRADIENT SEARCH ROOT IKE. 8E f SEARCH PARAMETERS i 
DKKAIILI VALUES ARE 1 

« UK STEPS. N- A STEP S1?K«D<* .iOOOOvK-Ol DOTH* . IOOOOOE-06 

wish fu changey 
r <> 

HEW LUST* .. I<<S2K2K4<<H 

COST FUNCTION’’ ..<«S2!5:. , K4«H 

wish io continue ihk starchy 

f A 

GRADIENT SEARCH RIIIIT I HP . ST f SEARCH PaRaMCTERK. 

DEFAULT Values area 

• OK SfKPS.N* » STEP SIZE.D* • iOOOOOK-Oi Dnlrl" • A<IOO<IOK <IA 

HASH TO CHANGE/ 

* l 


ENTER NKU VALUES I 
Mil >i . AE-S 


OKU COST* 
LAST SILP 
STEP SIZE 
MEW COST--. 
LAST STEP 
STEP SIZE 
NKU COST* 
LAST STEP 
STEP Sl/K 
NEU COST-. 
NKU COST* 
2 STEPS 
LAST STEP 
NEU COST* 
LAST SIEP 
STEP SI/E 
NEW COST » 
NEU COST* 
A STEPS 
LAST STEP 
NEU COST* 
NEU COST* 
2 STET’S 
LAST STEP 
NEU COST* 


.i2tt«v.sctov 

NOT ACCEPTED I 
REDUCED TO! . KOOOOOK-O I 
■ 44MSC45K40H 
NUT ACCEPTED.! 

REDUCED TO! .2SOOOOK-OI 

,.<1/20SE40H 
NIIT ACCEPTED I 
REDUCED TUI . »2!J«00E -01 

.;:vih:i.K4oh 

..<J/JOSET«H 

UITH PRESENT ORADIKNT AND IJMtN* 
NUT ACCEPTED I 
• 2V2UHSK40H 
NUT ACCEPTEDJ 

REDUCED TUI . A/3000E-02 

•2H/S/1E40H 
. 2V2HH!5K40H 

UITH PRESENT GRADIENT AND DM IN- 
NO T ACCEPTED.! 

.2H2SVIK40H 

.2HS/..‘iE40H 

UITH PRESENT GRADIENT AND DMIN* 
NOT ACCEPIEDI 
• 2//42/K40H 


NEU COST* 
2 SIEI'S 
LAST STEP 
NEU COST* 
NEU COST* 
2 STEPS 
LAST STEP 
NEU COST* 


.2/VKO/K40H 

UITH PRESENT GRADIENT AND DMIN* 
NOT ACCEPTED! 

• 2/2/1 /ETON 
.2/S.JASE40H 

UTTN PRESENT ORADIENT AND DMIN* 
H OT ACCEPTED! 

• 2A/AAHK400 


NEU COST* .22«.UIE4«R 

2 STEFS UITH PRESENT ORADIENT AND DOTH* 
LAST STEP HOT ACCEPTED/ 

NEU COST* . 2A.I24IE4 0H 

NKU COST* .2AA4.IAE40H 

2 STEPS UITH PRESENT ORADIENT AND DMIN' 
LAST STEP NOT ACCEPTED! 

NEU COST* .2K0AI2K40H 

NEU CCS I- .2AUIV0E40H 

2 STEPS UITH PRESENT ORADIENT AND DMIN* 
LAST STEP NUT ACCLPTEDJ 
NEU LOS I* .2S 44/ VET OH 

NEU COST* ...'K022.JK400 


:• STEPS WITH PRESENT ORADIENT AND DMIN. 
LAST STEP NOT ACCEPTED J 


. I.JSOOOE-Ol UERE I AKI N 

. A2ST<<<<<E -02WERE I AKEN 
.A2SO0OE-O2UKRE I AKEN 
< A2K000K- OIlUERE 1AKEN 
.A2S<i<K>E-«2UEIiE TAKEN 
. A2S<I0<IK <<2UERE TAKEN 
• A2S<I<I<<E 02UEKE TAKEN 
. A2"00«E-02UERE TAKEN 
.A2SOOOK 02UERL TAKEN 



j 


I 


! 

t 



\ 

| 

| 

» 

i 


I 


,4 


n 



A&V.W .lb 


t 
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Ntu tour- .ssosxvfiow 
cos; runt; r/ on* .asoaa/ooii 

wish io ooMrxtJtiK rut: search? 

■THE SEARCH IS CONTINUES ANS A riNAL COST FUNCTION IS CALCULATES.^ 

GKADXKMf SEARCH RIHJf IKE. SK f SEARCH EAftAHE TKK8 » 


CUAMLf VALUES ARE I 

« (IE (ifKI'S.M- 10 SI O' S1/.E.IJ- . A2S000E -02 


COST KUNf.TJOM* 
HA IRIX V ( 


.'fBSi 34E SOT 


. 40244531 1V023SE+0? 

.i:K!a53472i4V3?fe+<»i 
. 20.17.104144971 4E+00 

. 1244170:14439 ObK+02 
. j 27431 A»7SH104.»:*00 


4 .7978011 SO 10709ES01 

-.4V307i8014iVlM100 
HA IRIX VKS 

1 


,XV4S4iFI02 
. 67*321 OOi 

• ASivssr.too 


2 

.3 14 34 (JO 00 
.72S372E40I 
..13.:A3f;i:+O0 


2 

. 786744087272931. 10 A 
. t9»02U03020206Et02 
.3D436S»32tl042Of:iOl 
.54“29534274572*L'-0i 
3 


limit- . umowr-o:; 

3 

-.4141 VtlAO JVV V74t" +00 
-. 70040/344S1.13.1E.+00 
. iikciaooi .i/:m:s;»i +oi 

. 13998 J4S304AS4I: 401 


4 -.A/23S9E -01 .A/SOASirtOO 

wan; ro continue search? 

; « 

CAIN tlATRXX 

1 

4 


-. 1029031.100 
.7821821: Oi 
• aoiviioc+oa 

-.fm*79Eioi 


-.34024ii:-0t 
•9/40U4K 01 
. 1 13366L+00 
-•.lVl.17Al.IOt. 


1 . 131C0024UV4V4AE+O1 

-.»72.1A*.;«S7i3.l/Et01 

.! -.3GS3/9UU320347+.+01 

«3l 10401V7380* VI-102 
HA I R ( X AHrtfi 

1 


- . 14804 40 1 491 1.0 r. I(H 
. S033??»97All42?i: I oo 


I -.21 7/57+: +01 
- . 7423H4t:+00 
.1 .7435SV+.+00 

4 . :*32 v;*i.*e+oo 

RAUSE 1 :• 


r xKORHcru 


.3A424«E+0i 

.2V4204K-O2 

• isaVnE+oi 

-.11V43VEIOO 


2 

.. 1984001" +01 
-.1034.12+. 101 
.2HVA43I+00 
..ISAIOVL+OO 


2 

7 

-.31.1020+:iOl 

. 000000 ". too 

- .403424C I 01 
.OOOOOOEIOO 


.1 

-.1800421:+ 00 
-.314X721:100 
.M397V(jntOO 
.841H28C+00 

3 

- . 1 75A30+: f 00 
. 24V7U4I.+01 


. iAV.{?if*4o«:'4in:»«o 

- . so ion 72ii;'.:;v;>or:4o » 


. 2.100 721:101 

. 1201 l/TC+Ol 
-.I14977SCI01 
-.430791001 


•.371042000 
• 4.J4400M01 


-.137192000 
- . 773470+;- 02 



rwuil.' 
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MlsiTTU'ivrs^Qv t It.ij 

f « 

NS* / NO* 4 
-ALAN12TI 

i 


• M4000K400 - , A .I/000K400 


.200000E-01 


.OOOOOOKEOO 



2 ..I20000E-01 . U2000E401 .OOOOOOLEOO . OOO< 1 O< 1 K< 00 

* 

A . «S4<KKiE+<>0 -.24V000K+00 .VVAOOOEKOO ,0000001: 100 i 

\rou most select uirrcn j outputs u tu. m: ot;e» to run ntr iwm; : nvut f 

vie outputs corkespgnu to columns in -ai.ami 2 i 

SF.LKCIK.li COLUMNS UU.I. KOKH A MATRIX THAI MUST EiC NONCINGUl.AR 1 

ENTER OU'I PUTS TO BE UCED UNIF.OKRI ( 

» » ' 



tttiiitit*** i»7TiTfr7>> ; s X ON EiGETiVAl.llHT KOI; TJTj?Tt.RVi:i!****4***«* 

«*«« **************** mode 2 iki genual ur assignment *** ***************** 

ENTER Of; CHANGE riGCHVALUCSi 



BREW tons VhI lit:*, f 
» 0 

LAMBDA i I 

REAL* ISOOUOK.KOl 1KAII* . iftOOOOEIOi 

t , 

wish rn change* f 

• 1 . ! 

ENTER NEW VAt.MK(K) I ■ ■ • ■ '( 

* -ft O 1 • 

ML XT K TOT.HVALUK i { 

PREVIOUS VALIIKSV 

» 4 

LAMBDA 21 . 

REAl * •. IftOOOOETOl IMAfl* ISOOOOl 401 I 


WISH 10 CHANGE* . 

* I . . ! 

ENTER NEW VALUE* K) I I 

*-40. j 

NLX1 EIOENVALOEI 

PREVIOUS VALUES* > 

» O 

LAMBDA .11 • i 

REAL* -.2000004401 I MAG* .1000001401 

WISH TO CHANGE T J 

T 1 ! 

CNTt.R NEW VALUE* H) I 
t -7 0 

WISH 10 EXU ERIIM THTH MODE* 

» « . 

*************************** EXITING MODE 2 **»***«**»*«*> 



miiUMtllltT ASSIGN f. IGT.NVEC TGRfi ET/R ORKCRVCR************** 
**t*M******«*tM**t MODE .1 1 K11JTN VECTOR ASSIGNMENT *♦***»♦***«*«**: 
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ENTER OR CHANGE EIGENVECTORS l 


PREVIOUS VALUES? 

T 0 

EU.ENVKr.HIR V U <RKAt. ) 

.4H7000K402 
-•U/000K40* 
...... .1240004:402 

^UlbH TO CHANHKT 

ENTER A NEW DESIRED VECTOR J 
1 l 0 o 0 O 0 

deputed vectors 


(IHASO 
.AVA0OOF.4OS 
• SV1000E40;> 

hmvoookto* 


. A00<I<I<IK40A 

actual vector i 

. A00OO0K4O* 

e:rror vector i 


.OOOOOOK Too . OddOdOi; too 
• OO< 1 O<I< 1 E Kid . <i<i<i(i<aiE | <id 


. 000<IOOE4<1<I • OOOOOOK f dO 

I.E NOTH OK THE DESIRED VI CTOR . 

1.1 NO TH or THE PROJECTED VI CTOR:' 
LENGTH (IT fill: ERROR VECTOR „ 
^18 fHE ERROR ACCEPTABLE? 

NEXT EUiEHVECrOlil 
EIGENVECTOR V 2! < REot. ) 

.4II/000E402 

- • A.I/OOOKKH 


_ .4240001:402 

WISH 10 CHANGE? 

r t 

CNfl.R A NEU DKS/REU VECTOR I 
T O O A 0 O 0 

DESIRED VECTORS 


. 0<IO<I<IOEIOU 
I .OOOOOO 
I • OOOOOO 
'OOOOOO 


UMAO) 

* • AVOOOoEfOl 
-•AV4000i:i02 


IVOOoKIOi 


< ooooooe too 
ACTIIAl. VECTOR | 

•OOOOOOE40O 
ERROR VECTOR I 


• KlOOOOETOl .OOOOOOK TOO 

• I OOOOOETOl .OOOOOOK 1 00 


• 000<I<I<IE400 . 0O<IOO<IE4<IO 

length of the desired vector » 
!'K N 5™ ,,K f,,E E'iOJtC TEN VECTOR- 
LENGTH or THE ERROR VECTOR « 

S8 THE ERROR ACCEPT ARI.K? 
f i 

NEXT EIGENVECTOR! 

EIllENVECfOR V .Tl *<REA1> 

-..I04000E400 
-.VOS 0<I<IE40<I 
....... .479000K40A 

WISH TO CHANGE? 

? S 

ENTER A NEU DESIRED VECTOR I 

f 0 O 0 <l ( <i 

desired veciors 


. <I<I<IOOOE 100 

1 . 000000 

5 . 000000 
•OOOOOO 


< THAO) 
.J44000K400 
-.V70000K-01 
.70VOOOK-OA 


. 00<I00<IE4<I0 
ACTUAL VECTOR I 

. <I<IOOOOE4<IO 
ERROR VECTOR I 


. 0<I<I00<IE400 • 100000E40I 

.OOOOOOK 400 . 100000E40S 


«•« ittpzs*, tssnr -«sr 

888! s 85 SKIIZ"": ‘:«rs 

SH THE ERROR ACCEPTOR! E? 'OOOOOO 

t t 

^R'xTr ,,N,EMfS W ' ,:UUkNt ' W *TA EIRE INCLIITIKI 

l 


.sooooooooo«oooE4<u .ooooooooooooooetoo .oooooooooooo<ioeioo 


I 


OtiiUlNPL g 

Of V 1 0 0 i'i QJAw»< ■ 


. 0<JOO<JO<JOOOOOOOK400 


• ioooooooooooooi: jox 


i • <J<J<J<JO<J<JO<J<J<JO<J<JK4<J<J , <JO<J<JO<JO<J<J<JO<J<J<JK fo<J 

UISH TO IHSRI.AY THE NORMALIZE)) E T OKNVKC TORS 7 
i « 

<)A<N MATRIX Kf 

I 2 


1 . 1V018A2S27/VAVE4O2 

s -<saa24443as433vk4oo 

.< -<320sbu72043ia24moi 

HAIR IX AHA T 1 

1 

J - . soooo(>o<»<k><>{>o t r + <• j 

2 <J!iS:?ynA7OH0OB4>14 


. 7AABSA1 X 342V WE 1 0 X 
.1A7452AVA01H7/1402 
. sooioooobabbosi'. to x 

-..UV74423X0V205L 13 

- . A<I00<HI<I<I<100002IM 0 X 
<2U42i70V430404i>i3 


< ooooooooooooooe4oo 

• 100000000000004:401 


- . 10B71 A/<J27704AE400 

<?>'.50aivxsoaoiovk: 02 

-.AVVH0/0AAB7A07K40J 

.1 

. X 77A3%.V43V 400.SC - X 4 


-.SMSmSJ2;U.*»M.XA 
•• • /oooooooooooooi: i ox 


W1BH to «H?”ho!S 4 ?M?5 # i5Sw? • 20421 70V430404C- 13 -./OOOOOOOOO 

T 1 • 

4****t*******M*ttt**ttttt* eyjTjm wmt j ************************* 


i; 

LI 

1 

• 

2 

3 

4 

i- 

f 

i . ivute<se4o;> 

-<3AA24AE400 

<32<JBBVE4<J1 

! 

<<J<J<IOOOE 400 ' 

<“ 

7 ./AABSAlMOt 

< 1A74BBE402 

. 3<I<H<1HIM<IX 

, 0<J<JOO<JK 400 

* 

.1 1BB/174MOO 

.S30A1VK-02 

<AVV<J<J7E 401 

<O000<J0E4O<J 


MATRIX Ei 

X 

2 

3 


X B00<I<1<IK4<U 

..IBB27XE-14 

.X42I0VE-13 

• 

4 

2 -.31V744K-X3 

A000<J<JE4<J1 

.2B4217E-13 

J 


J . 1/7A3AK-14 

-.BS5112E-1A - 

.7000004; 101 



MATRIX kt 

1 

2 

3 

4 


X .2BAA2V4MOO 

-.2V2BV4K-01 

• 203 USE 402 

-.371BA2E400 

i 

2 < 1843/HE 400 

.71 0B7BE400 

.1713H7K402 

< 

<434AB<JEI<J1 

j 

PAUSE 20 

r 

MATRIX (11 

1 

2 

3 

: i 

! 


j -.»yx«;»HE402 

.274 7204.401 

• X7B10BiM<J2 

i 

.12374 74,400 | 

5 

2 -.421A0AE402 

-.b;tavb;»E402 

.57A1SJE407 

1 

-.X1SB42E400 


.4 -.B2V2H4E401 

-<34ahiac-oi 

.4B4B13LM02 

<77<J12AE40<J 

ct • 

MATRIX TBI 

1 

J. 


/ \ 

X . 2<JO<I<J0E4O2 <0000001:400 

2 <00<J000E40<J . 250000E402 

.5 , <JOOO<J<JE 400 <0000<J0E4<J0 

E4I TER _ IFiTTXAL CONIU TIONS FOIT OR «i <NAI 

nTTATEC 



r o o o l o o o 
USE THE KOI.I.OUIHO 

INITIAL CONDITIONS 

IN 1IM4: REM'. 



. ooooooooooooooK4oo 
. OOOO<1<IOOO<1OOOOK4OO 
.UOOOOOOOOOOOOOE 400 
> XOOOOOOOOOOOOOE40I 


I 


1 


ORIGINAL PAS 2 H *3 
OF POOR QUALS 7 V 


S . OOOOOOOOOOOOOOE+OO 

* .OOOOOOOOOOOOOOETOO 

'/ ,OOOOOOOOOOOOOOK+O0 

H , OO< 1 < 1 O< 1 OOOOOOOOLKOO 

f . OOOOO< 1 OOO< 1 <I<K 1 OKK OO 

10 .OOOOOOOOOOOOOOEKOO 

f M18H TO DO f XML' RESPONSE* 

«Mt«****tt«*«4****t*t* H(m - „ n , ML . KJ(m(1AU0N 

********** CHOUSE SIMULATION OPTIONS! 

^ -ENTER! X TO SIMULATE KA.T, 3 TO SIMULATE LAHAf.I.CJ nil} I.AlU.I)! 
^ENTCK 0 TO SIMULATE OUTPUTS. 1 10 SIMULATE STATE VARIABLE!.! 
^ENTLR SIMULATION TIME > <KEAL. NUMBER IN SECONDS >i 

«™ N,Wfl,;R " K R0 ^ f ' ; ro be calculated, < aoo maxm 


SPICJFT THE INITIAL CONDITION!;! 

X 1(0)1 
f 0 

X ;:<o>i 
t 0 

x :»(«)! 

TO 

X 4(0) 1 
T 1 

X S(0)l 
T 0 

x 6(o>: 

T 0 

x /((»: 

TO 

x H(o>: 

T 0 

x v( 0 ): 

T 0 

xio(o>: 

TO 

CHOOSE INPUT OPTIONS) 1 EUR NO INPUT. 3 EOS A STEP INPUT. 

INPUT OPTION EON U l.’’ ^ 4 K0K A TR " Mf;A,KD «*«>'< 

T 1 

INPUT OPTION PUR U 21 
T 1 

^ENTEN 0 FOR 80 DISPLAY COLUMNS. 1 FOR 12V COLUMNS: 

^ENTLR 0 FOR INDIVIDUAL AND 1 FOR MULTIPLE PLOTS! 

^DO YOU WISH TO SEI THE MIN-MAX KANSK S FOR THE AXEST 

ENTER MIN X.MAX X.MIN Y » AND MAX Y VALUES! 

*0 12 -*00S *0415 

POSITION PAPER AT TOP OF FORM AND TYPE ANY INTEGER 
YOU MAY ADD A SHUNT NOTE (20 CHARACTERS.) 

THE TIME RESPONSES ARE NON PLOTTED* RESULTS ARE SHOUN IN FIGURES S-lt-S-ES* 


/ 








m 




»%; r -3 ;,'.', ; V,:'; > , 

■■jPIft 

iBfci* »v Ti *Xtf‘i'$-.i'f i: V> -'". 
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